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Expansion Procedures and Similarity Laws for Transonic Flow 


Part I. Slender Bodies at Zero Incidence 


By Juxian D. Corr and Artuur F. MessiteEr, Pasadena, California!) 
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1. Introduction?) 


The purpose of this report is to provide a detailed and comprehensive 
account of a transonic approximation as applied to flows past wings and 
bodies. It is mainly concerned with the derivation of approximate equations, 
boundary conditions, etc., rather than with the more difficult problem of the 
solution of transonic flow problems. Thus the report contains for the most part 
a re-examination of the basic ideas, as presented for example in [1]). The 
essential new point of view introduced here is to regard the approximate 
transonic equations as part of a systematic expansion procedure. Thus, it 
becomes possible, in principle, to compute the higher terms of this approxi- 
mation or at least to estimate errors. 

In the next section the form of the expansion and the reasons for it are 
explained. In the succeeding sections the equations of motion, shock relations, 
and boundary conditions for the flow problem are presented and then the 
expansion procedure is applied systematically. ' 

The resulting system of equations for the first, second, and higher approxi- 
- mations is presented in section 5. The main results of interest for practical 
applications concern similarity laws and the pressure coefficient on the surface 


1) Guggenheim Aeronautical Laboratory, California Institute of Technology. 
2) This research was carried out under ARDC Contract AF- 18(600)- -383. 
3) Numbers in brackets refer to References, page 29. 
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of slender bodies and these appear in section 6. The remaining section treats 
bodies of non-circular cross-section. 

The authors wish to acknowledge the assistance of S. SHEN in the prepa- 
ration of this report. 


2. Motivation 


Linearized theory (for a complete account see the book by Warp [2)) is 
useful for estimating the pressure on slender bodies and wings in the subsonic 
and supersonic regimes. However, the results of the theory [2] predict an infinite 
pressure for bodies flying at the sonic velocity, Mach number M,, = 1. This 
result, so far removed from reality, indicates a failure in basic assumptions 
underlying the linearized theory. It also indicates that the linearized theory is 
a poor approximation for some range of Mach numbers about 1, the transonic 
regime. It is the aim of transonic theory to provide the proper mathematical 
description of the flow in this range and hence to allow calculation of pressures 
and over-all forces. 

From the point of view of mathematical expansion procedures the ‘break- 
down’ of linearized theory corresponds to a non-uniformity in the expansion 
procedure. The successive terms in the expansion, which are implicitly assumed 
to be small, become comparable with the first term (which is given by linearized 
theory). 

In more detail, in a typical case there exists a velocity potential ® which 
for steady flow past a given body depends on the dimensionless space coor- 
dinates (x, v) and the physical parameters 6, M,, y: 

6 thickness ratio= 7/Z , 

# maximum thickness of body , 

L characteristic length, e. g., over-all length of body , 
M., free stream Mach number = U/a, , 

U flight speed , 

a. speed of sound at infinity, 

y ratio of specific heats. 


@ is determined as satisfying a non-linear partial differential equation in 
(x, 7) embodying conservation of mass, momentum,.and energy, and certain 
boundary conditions. The body causing the disturbance to uniform flow is 
assumed to be thin so that the potential may be represented by the following 
type of expansion 


P(x, 7; Moo, 0) = UL {x + e(3) g(x, 75 B) + €0(6) galx, 7; B) +--+} (2.1) 


where B =V| 1 — M%,|. The ¢,(6) > 0 as 6 + 0 and e,,,/e, > 0 as 6 > 0. The 
first term in the expansion represents a uniform flow toward the body and the 
second term ¢, y, represents the potential of linearized theory. The remaining 
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terms & Po, €3 PY; etc. represent second-order-and third-order corrections etc. 
In any given case the e,; are known functions of 6-determined from the equation 
and boundary conditions. By requiring the non-linear equation for ® to be 
satisfied identically for all orders of magnitude in 6, the linearized equation for 
@, is found. Typically it is of the form 


(L— Me) Meet Mp, te = 0. (2.2) 


In order for the expansion (2.1) to be a good approximation it is implicit that 
the velocity perturbations due to successive terms in the expansion are of 
successively smaller size over most of the flow field. Unfortunately this is not 
possible when 8 > 0, M,, > 1 because P1, > ce; thus the expansion is said to 


be non-uniform with respect to 8. That is, for any given , a sufficiently small 
N(f) can be found so that for 6 < N(f) the approximation is good‘). But as, 
Bb +0, N +0 also. 

That such a non-uniformity might exist is evident from the mathematical 
structure of (2.2). For M., <1, , satisfies an elliptic differential equation simi- 
lar to LAPLACE’s equation, and for M,,> 1, q, satisfies a hyperbolic equation 
similar to the wave equation. Due to the completely different dependence of 
solutions to these equations on the boundary conditions a non-uniformity 
might be expected. It is well known [3] that according to the exact equations 
a flow which is locally supersonic is described by a hyperbolic equation (has a 
wave structure) and flow which is locally subsonic is described by an elliptic 
equation. In the transonic range the velocity may change from subsonic to 
supersonic in the flow field. Linearized theory is too, crude to describe this 
feature of the flow but one will require it of a transonic theory. This implies 
that the basic equation of the transonic approximation is one of changing type 
(elliptic to hyperbolic). This also implies that the equation is non-linear as the 
line across which the change of type takes place, the sonic line, may not be 
fixed in advance. Thus the transonic expansion must have a substantially differ- 
ent form from the ‘linearized’ expansion (2.1) and in particular must involve 
different combinations of parameters M,, 6. 

Before outlining the proposed expansion some physical aspects of the 
approximation will be discussed. Linearized theory is identical in content with 
acoustics; squares of disturbances are neglected (for a complete discussion 
see [2]). In a system of coordinates fixed in the air at rest at infinity the flow 
is described as a system of acoustic waves produced by the motion of the body. 
These acoustic waves travel with a constant speed a, and each wave carries 

a small disturbance. When a body flies with a speed U close to a,, some of the 
"waves emitted by the body are always close to the body and a build-up of 


4) Of course, good has to be defined more precisely. For example, it might be demanded that 
the first omitted term ¢, p, amount to less than 10% of & 9. 
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pressure on the body may result. This result depends on the linearized descrip- 
tion of the flow and causes the failure of approximations near M,,= U/a,.=1. 
Actually, as these waves build up their interaction becomes important. The 
result is the famous steepening of the wave front described by RIEMANN ((3], 
p- 45) in his theory of plane waves. In RieMann’s theory it is shown that non- 
linear terms like u u,, are responsible for the steepening of a plane wave travel- 
ing in the x-direction (u = flow velocity component in x-direction). From this 
point of view such a non-linearity is to be required in the basic equation of the 
transonic approximations. It is, in fact, the same non-linearity discussed in the 
last paragraph. 

Another physically significant point concerns the relative orders of magni- 
tude or perturbations from the uniform stream. It is implicit in the linearized 
theory that the (dimensionless) velocity perturbations in the x- and 7-directions 
€1 Pi,» &1 Yi, are of the same order of magnitude, namely ¢,. On the other hand 


an examination of the exact solution of a simple supersonic flow [3], the 
expansion around a corner, shows that the perturbations in the x- and r-direc- 
tions (u,v) from a uniform sonic stream in the x-direction are related as 
v ~ uv”, A similar result holds for shock waves in the transonic range. This 
sort of order of magnitude relation will be required, in the large, for the per- 
turbations of the basic transonic approximation. If it is known that a potential 
exists (as will later be shown to be the case) this last requirement has the 
interesting consequence that the expansion must be carried out in a distorted 
set of coordinates. For example, represent the potential by ¢,(6) (x, 7) where 
v = eso that the actual velocity perturbations 0/0x (e, 9), 0/07 (e, y,) are 
of the correct order of magnitude. A more vague interpretation of this coordi- 
nate 7 is given by saying that in the transonic range perturbations caused by 
a body in the flow extend much farther in the direction transverse to the flow 
(large v) than in the flow direction and that a shrunken coordinate 7 must be 
used to bring these disturbances back in view. 

All of these physical considerations are keys to the more general form of 
expansion which must be used in the transonic case. The idea of the expansion 
is again that small disturbances in a uniform stream are caused by the intro- 
duction of a body of small thickness ratio 6 and the expansion procedure is 
based on 6 > 0. However now a more general form of expansion, allowing for - 
different parameters and distorted coordinates must be used. Assuming a 
potential, a sufficiently general form is 


O(x,7; M,., 6, y) = UL {x + 6,(8) p(x, 7; P, Pp) + €9(6) po(x, 7; PB, P,) +-- 2 
(2.3) 


where P,, P, are functions of (y, M,,, 6) which are independent of 6and 7 = Ar 
where A = A(d, M,,, y). 
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The expansion parameters ¢,, P,, P, and the coordinate distortion A must 
be determined in the course of finding the equations for y;. The main require- 
ment to be imposed is that @, satisfy a typical transonic equation. The choice of 
parameters in such an expansion is of course not unique but it will be shown 
that all such transonic expansions retain some typical features in common. 
In the next section the basic equations of motion are presented and in the 


sections following that a general type of expansion procedure corresponding to 
(2.3) is applied. 


3. Basic Equations of Motion and Boundary Conditions 


The basic equations of motion consist in statements of the conservation 
of mass, momentum, and energy in the form of differential equations and jump 
conditions across shock waves. To these must be added an equation of state 
and the Second Law of Thermodynamics. The medium is assumed to be a 
perfect, inviscid gas. The region under consideration always extends to infinity. 
The boundary conditions prescribe the flow at infinity and require the flow to 
be tangent to the body surface. 

Assume that a characteristic length L exists, say, the body length and make 
the space coordinates dimensionless with L. Then the differential equations of 
motion can be written in terms of dimensionless coordinates as: 


continuity divog=0, (3.1a) 
Ss gq? Eee 1 

Euler q-Ve= v(4) +0-G=——P, (3.1b) 
= P 

entropy q: V(y) 708 (3.1¢) 


System (3.1) together with the shock relations, which will not be written out 
((3], p. 51) comprise the basic system. An expansion procedure could be applied 
to the basic systems; however here a modified system will be derived which 
contains some integrals of the basic system. 

First: note that (3.1c) implies that P/o” is constant along streamlines, 
although of course the constant jumps across a shock wave. For all problems 
considered here, there will be a uniform state at infinity, so that (3.1c) implies 


=e =k (say) is constant along a streamline. (3.2) 

Q” C3 

The basic integral invariant is formed by taking the inner product of g and 
the Euler equation and integrating along a streamline. The result is 


q? a Peas cage eeamnes | A Guan 
: +/ See agen const along streamlines (3.3) 
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Furthermore, the shock wave relations show that the constant in (3.3) is 
conserved across shock waves. Since a uniform state is assumed at upstream 
infinity the constant for each streamline is the same. Thus, the following 
invariant is obtained for the entire flow field 
a? U2 ies i (24. 7) ak? (3.4) 


q* uae be Past as 
Drie esl tae i | 2\y—1 


By forming the gradient of the invariant (3.4) the vorticity laws are derived. 
Replacing V(q?) from the Euler equation (3.1b) and using (3.2) gives the vortex 


law. 
= : VR 345 
OG == s ( . ) 


Finally, the density may be eliminated from the continuity equation (3.1a) 
by evaluating g x V using (3.2) and (3.1b). The result is 


Pde gs a: v(F) (3.6) 


Summarizing, the original system of differential equations can be replaced 
by (3.4), (3.5), and (3.6). 


continuity a?divg=q: o(2) 3 (3.6) 

: 2 q? a? eS ; as 1 /y+1 

invariant => + ae a ee ee ae 7) ae?) (ana) (3.7) 
Se Se Ss eo” . 

vorticity o-¢g= at Vk, (323) 

and q-Vk = 0. 


In certain applications, for example in flow which is subsonic everywhere, 
the entropy and # are constant throughout. The vorticity law (3.5) then implies 
@® = 0 and a velocity potential exists. Then (3.6) becomes the single equation 
for the potential. It will be shown later that for the transonic flows considered 
this statement is true (to a certain order) even if shock waves are present. 

It is convenient later to use certain relations deduced from the basic shock 
relation. These are the shock polar equation which relates the velocity compo- 
nents, and equations for the pressure jump, wave angle, and density jump. 

The shock polar can be written 


2 PRO KC 


Polar vy)" = (y — u)) (3. 8a) 


2 2 2 
(0) Soe WO) 
et u ur’ “ud! — a* 
where superscript (0) denotes conditions ‘upstream of the shock wave and 
superscript (1) denotes conditions downstream. 4 is the magnitude of the. 
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velocity upstream. x2) ss the component ofthe velocity downstream in the 


direction of «') and v™ is the component normal. u,v are related to the 
Cartesian components as follows 


0 2 2 
yO 4? + go | 
(0) (0) 
wf) = g@) ___ Lee Dk : ath | 
V2 + V. OF 4 gi) (3. 8b) 
0) (0) 
ee ee, de heieesiss | 


qx aie a ae = 02 q; cen me Vc . 
Va? + a® Vasa? | 


A special case of interest occurs when the shock is in the free stream so that 


wu — U and is in the x-direction. For this case the following relations hold 
0 
Wave angle tan 0, = —=——__—_, (3.9) 
berate 
pressure Pe = Pee UU) ge). (3. 10) 
(1 (1) 
density a) i ten 6. (SL) 


In case the flow does not have axial-symmetry these formulas still hold if 
qg, is replaced by the resultant transverse component behind the shock. 

The boundary condition to be applied at the body surface is one of tangent 
flow at the surface 


qg:Vs=0 where s(x,7,0)=0=7-— 6F(x, 6) (3.12) 


is the equation of the body surface. Written out (3.12) is 
q,(*, OF, 0)  qy(x, OF, 0) , OF _ qo(x, OF, 0) (1\ (OF 
ee, SaaS i (=) (oo) 


(3. 13a) 


4 Ox U F} \ 00 


and for the special case of axial symmetry 


oa = pe 6F’(x) where on the body ¢ = dF {a (3.13) 


The other boundary condition is a statement that all flow disturbances die 
out at upstream infinity. 


4. Deduction of the Form of the Expansion; Body of Revolution 


In this section it is shown how with the aid of certain general requirements 
the form of a transonic expansion can be deduced. The requirements are stated 
for a free stream Mach number of one (M,, = 1). These requirements are: 
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(i) The shock wave relations should not degenerate in the expansion. 
Rather the shock waves should approximate real ones with jumps of velocity 
and pressure etc. 

(ii) The differential equations of motion should be able to describe the 
flow satisfying boundary conditions on the body and at infinity. These two 
requirements define the expansion in terms of 6 uniquely at M,, = 1. However 
the extension of the expansion for M,, + 1 is not unique, in a way which is 
shown in detail below. It is important to note that these requirements yield a 
transonic approximation which has all the qualitative features described in 
section 2. 

The shape of the body of revolution and the boundary condition on it are 
given by equation (3.13b). The expansion procedure is based on 6 -> 0 and it is 
assumed here that (M,, 6) are the parameters defining the flow. In accordance 
with the remarks made in section 2 the following form of the expansion is 
assumed for the velocity components 


i, %; Moo, 6 A . 
qa ¥ oy eet ++ &(0) u(x, 7; Py) + &0(0) u(x, 7; P,)) +---, 
( Mp, 8) (4. 1a) 
why, V5 oo? is 2 
aol 7 = v,(0) v4(x, 7; Py) + r9(8) vox, 73) to, | 
where the e; and y; each form a decreasing sequence as 6 0 and 
y= 6%r, P, = parameter independent of 6. (4.1b) 


In accordance with requirement (¢) the assumed form of the expansion (4.1a) 
is substituted into the shock polar relation as given by (3.8). Thus 


(0) 


U 
a Elta ee, 7 = lei aan ae 
(4. 2a) 
on (1) (0) 
U = 9, (vp — oy’) + , 
and the shock polar, for M,, = 1, a* = U, becomes 
A 2 (0) (1) ess =— 
om (vi) — yO)? 4... = 2? (yo — uid”) ae & (wy) + uy!) + 2 
[iS ES eee 
JET ee hr eel 
or 
2 /.,(1 1 
vf (of — of)? + veo ob (al! — ol 8 fa + oA ZO 4 2b) 


In order to avoid a degenerate polar it is necessary that yi~e? as 6>0. 
Hence we take ; 


v= eh. (4.3) 
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The assumed form of the expansion is now used in the modified continuity 
equation (3.6) and use is made of the invariant (3.4). Thus 


q? 


ie =14+24%44+28u,+ 424 2&4 Ey Uy Wy + ef UZ +++. + ev? +... 


’ 


a® 1 et : (1 ce 
4 5 4 


= ae ~ Y — Vey — (y — 1) gw, — 2 a 


7q 0 3/9 0 
Via (Lt ey Wy + by ty + --) Bee Or ee 


and equation (3.6) becomes 


it SN Sex ous 
Lae — 0 — 1) ey — (y — 1) ep ty — 2S * tag —--] 


19 1 ib 


= [(L + ey + ep tty toe) sei? Ho ge te] i 


2 2 
2 “4 EN Gh ae 
x [en a + 0 te + of E+ oy ep ty ty + uh + 
ef UY 3/2 


neato tee My 0% +o]. 


The form of equation (4.4) when M,, = 1 is now studied: 


[2 = =V) em — (y— 1) eee ys he | ] 
il 
x [ex t, + 22, + °° + 8? 60 (4, + = 1) + v 6% (+ = jen] 


0 9 0 
_ [(1 + ey + 82 Me + oo) Sp Fat? Oty Ge +e (4.5) 


| 
uy 2 U3 
4 [ex tr + be Met Po + ey ety Me tee D to 
ef Ui 3/2 


ame eon et ¥5% % + -+-| 


and the requirement is made that the largest terms in equation (4.5) yield an 
equation which contains the possibility of describing the flow field in the large. 
It is easily seen that the largest terms are 


it 
esl? 4 (a, =i oe 04] =e(y+1) my My: (4.6) 
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If either 5¢ = o(et/”) or el? = 0(6%) this contains only one of the unknowns 
and yields unreasonable relations. Therefore choose 


ge age (4.7a) 
and obtain the first approximation equation for M,, = 1. 
eto n= +1 ym. : (4. 7b) 


The extension of this form of expansion to M,, + 1 can be performed in 
various ways by requiring the term e, [(1/M%,) — 1] #4, in equation (4.4) to be 


of order ¢7; that is 


il 2 
lee 1) &, = Olt) a8 seal (4.8) 


Note that if 
= o(| a2 - 1] é,] as & >0 

then a degenerate form of equation is obtained while if [(1/M3,) — 1] e, = o(e7) 
a more restricted approximation scheme is obtained. The requirement of 
equation (4.8) can be satisfied in a variety of ways, of which the following 
simple one will be used in this section 

1— M%, 

| 


= P,= similarity parameter independent of 6 = K (say). (4.9) 


As &,(6) > 0, M,.(6) is taken to be such a function that P, is independent of 
6; or 


M2,=1-—24(d)K. (4.10) 


The use of this parameter also extends the shock polar (4.2b) away from 


M., = 1. Thus, the largest terms in equation (4.4) yield the first approximation 
equation 


1 
Kw ++ > 4=vtiym,. (4.11) 


Next, another relation between (e,, 6) is found from the boundary condition 


of tangent flow equation (3.13b). Using the assumed form of expansion (4.1a), 
equation (3.13b) becomes 


eq!” 0, (x, 0°*? F) + vg 09(x, 6°41 F) + «+» =[1 + 8, u(x, 0°! F) + +++] OF'(x) . 
(4. 12) 


Equation (4.12), which holds on the boundary gives another expansion as 
6 ->0, the ‘expansion on the boundary’. If we assume a > —1, then the 
behavior of v,, v, etc. as 7 +0 is decisive in fixing the form of this expansion. 
If the physically reasonable assumption is made that u, is not too singular, 
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O(log) as y > 0, then from equation (4.6) itis seen that v, may be ~ 1/r. On 
physical grounds this type of behavior is to be-expected as it is recognized as 
the fluid source line on the axis often used to represent a body. Thus we write, 
near the body, 


Sie - 
Vy (x, 7) > el as 7-—->0O (S,= source strength) . (4.13) 


Our later considerations will show that the case a = —1 can not lead to any 
transonic approximation, while a < —1 is no good for obvious reasons. 


By using the limiting behavior expressed by equation (4.13) the expansion 
on the boundary (4.12) becomes 


19 Ont 
82 (2h a +) iS als FP (4.14) 


Equating the largest terms in (4.14) as 6 > 0 shows that 


a= oe (4.15) 
and 
Sy (ep (x) (2) (4.16) 


Equation (4.16) is the conventional result for source strength in slender body 
theory, namely the source strength varies as the rate of hones of the body 
cross-section area [4]. 

Thus the parameters for the first approximation are determined from 
equation (4.3), (4.7a), and (4.15) to be 


cet ey = 07 “9, 0" (4.17a) 
and 
1— M2, 
IO 35 (4.17b) 


K is practically the similarity parameter introduced by KARMAN in [1]. 

To complete the specification of the problem for the first approximation 
the vorticity laws which depend on the entropy gradient must be considered. 
An estimate of the entropy changes can be obtained from the shock relations. 
It is sufficient for our purposes to consider the shock in the undisturbed stream. 
Similar results will hold in the general case. The expansion (4.1a) now reads 

fe. =1+ 62 u,(x, 7; K) + (6) u(x, 7; K) + 
7 = Or (4.18) 
Ae = 88 v(x, 7; K) + 1918) vole, 75K) + 


and using this in the expression ee 10) for the pressure ratio across the shock 
wave shows pu 
Poo 


=1— @yu)+---. (4.19) 
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Similarly expansions for the wave angle and density ratio can be obtained 
from equations (3.9) and (3.11) respectively 


(1) 
tan Oy = — X Gar too (4. 20a) 


Consistent with (4.20a) the shock shape can be written 
Pi O4(a) a (4. 20b) 
the density ratio is 


ra La pee 62 wi) +... , (4.21) 


Superscript (1) means evaluation of the functions on the shock as given by 
(4.20b). Now 


oF Sal eye eee Ey ae 

n Gee (4.22) 
BR 
r= 1 + 0(6%). 


That is, at the shock there is no entropy change of order 6”, and since the 
entropy is constant on streamlines in the flow downstream of the shock a 
relation similar to (4.22) can be assumed to hold everywhere in the field. If 
there is more than one shock including some in a non-uniform stream the order 
relation specified by (4.22) still holds. Hence, in the entire flow 


Sa ae (4.23) 
Using this result in the vorticity law (3.5), the flow can be seen to be irrotational 
in the first approximation 


ae OG aunt 4 Ode 1 RO Ok 
U \U or foe) = ees ore 


— (1+ 624, + -++) [65(4,- — a ) —+-] 


LJ Gee hee 
vad (Leb 7 eee 


0 
é OF [o (0?) ] F 
so that 
This completes the specification of the problem for Uy, 0, and the first approxi- 
mation. 


To determine the further course of the expansion certain steps must be 
repeated for the next order terms. Returning to the modified continuity 


Vol. VIII, 1957 Expansion Procedures and Similarity Laws 13 


equation the second order of terms in (4.5) might be 


1 : 
Vo (v2, = vs] = (y + 1) 6? eg (1, ug), + 68 vy (uy. + %4,) 
(4.25) 
yp Lows ‘ i 
+ 66 Go yt My + 08 (y — 1) wy (4, apa v4] ‘ 


Reasoning as before equations capable of satisfying the boundary condition 
at the body and infinity will be found if 
Vo = €> fo) . (4. 26) 


a 


The extension of the equation for M,, +1 follows automatically now by 
MS, = 1 — K 6%. We have, as a possibility, 


2.0? K ty, + 28 (v9, + & vs] = (y + 1) 62 &y (uy Me), + O(8%) . (4.27) 


The second relation between e, and yg can be found from the expansion of 
the boundary condition on the boundary. To carry this out it is necessary to 
know the first few terms of the expansions for 1, v,. Introducing a potential 
Q(x, 7; K) to satisfy equation (4. 24) the equation (4.11) determines 9, 

L 
M+ S Pz-=[yt 1) m, — K] m,,- (4.28) 


~~ 
rr ae a 


The beginning of the expansion of g, near 7 = 0 and near the body has already 
been determined in equation (4.13). The first two terms 


Py = Sy(%) log 7 ++ 24(%) <e-*s 


of the expansion are solutions to equation (4.28) with the right hand side equal 
to zero. S,(x) is known but g,(x) is an unknown function of x. The successive 
terms in the expansion can be obtained by iteration using the first terms in 
the right hand side. The result is the following asymptotic expansion near 
the axis. 


= 1 , " Pe ~ 
@ye 7; = Sy(x) log7 + 2,(%; K)- (7? log?7 ) (7 _ Si St) + O(7? log7) | 
(4.29) 


0%, 7; K) = ais + O(7 log?7) . | 


Now two possibilities with respect to equation (4.27) need to be investi- 
gated, ¢, = 64 or 64 = o(e,). In the latter case the terms O(6°) drop out of 


equation (4.27) to give 
il 
us Us, ae Vp. ae FF Vo = (y + 1) (uy M2), - (4.30) 
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Then, applying the same reasoning as before 
V9(x, 7) > +--+ as 7 >0 (4.31) 


and y, can be found from the expansion of the boundary condition (4.12): 


43 [ 12) + 0(6%10g28)| + (S20 ++.) 


= (1+ 262 S/logd + ---) 6F’(x) . 


| (4.32) 


The second order terms in this expansion are thus 0(6° log 6) and the boundary 
condition can be satisfied if 


V_ = 6° logd (4. 33a) 
and, from (4.26), 
&, = 6* logo. (4.33b) 


This result is consistent with the assumption made that 64 = o(e,). If on the 
other hand the assumption is made that 


il a = 
So = 0", Vg ~ > logy as 7 >0 


instead of as shown in (4.31). This follows from the behavior as 7 > 0 of the 
terms of O(68) in (4.27). No consistent expansion can be found this way so 
that the possibility ¢, = 64 is ruled out. As the procedure continues the terms 
which were left out appear in the equation for (ws, v3). In fact S, can not be 
found until they are considered. 

The form of the expansion up to the second step is thus determined and 
the vorticity laws provide the second equation as before. The shock relations 
and the entropy jump across the shock must be considered first. A consistent 
way to treat the shock waves is by satisfying the shock conditions on the shock 
front 7 = o,() as found in the first approximation. In all transonic problems 
involving a shock wave the position of the shock cannot be fixed in advance 
and the determination of o,(x) is part of the solution of the first approximation. 
Some details will now be given. 

The wave angle equation (3.9) indicates 


faa 


tanO, = —-¢ X —G + O(6 log 6) (4.34) 


y) 


in view of the second terms in the expansion. This implies a shock shape of 
the form 


7 = 0,(x) + 6? logd o,(x) +--+. (4.35) 
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The velocity components behind the shock can be expanded as 


1) ‘ 
Ge 
Tr =o 1 + 6? u(x, o, + 6? logd a, + +++) + 6* log dé us, (%, Oy tees) eee 


(1) 


es < 

7 = 1 + 6? u(x, 03) + 64 log 6d [ u(x, 04) + 04 My~(%, 0;)| +++, (4. 36a) 
q,” 3 
a= 0° v4 (%4, 0;) + 6 logd [v,(x, 04) + 05 Vy~(X, 0)] +-°-. (4. 36b) 


Now, the pressure and density relations (3.10) and (3.11) become to the second 
order 


Pw : 
| xian 1 — 6? y u(x, 01) — 64 logd y [u(x, 01) + 05 Uy~ (x, 0z)| +++, (4.37) 
aor = 1+ 6? u,(x, o,) + 64 logd[ w(x, 04) + 02 Uy~(X, 04)] +++. (4.38) 


Thus the entropy change across the shock is specified by 


RO) 
Again, consideration of the vorticity law (3.5) shows that the flow is irrota- 
tional to the second order 
Up.— Vp =O. (4.40) 


The formula for the shock polar to the second order can also be worked out. 

The procedure just outlined is thus complete and can be carried out to any 
order. The results about the entropy jump across shock waves which were 
used in this section were derived for points off the boundary. Similar reasoning 
can be applied to shocks on the boundary to show that the flow is effectively 
isentropic there. 


5. Summary of Equations Derived in the Expansion 


Continuing the procedure started in the previous section leads to the 
following expansion 


fe <1 + 8? my (x, 7; K) + 41086 te (x, 7; K) | (5.1a) 
+ 64 us (x, 7; K) + O(68 log?6) , 
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fe = 6° v, (x, 7; K) + 6° logd 0, (x, 7; K) (5.1b) 
+ 65 vs (x, 7; K) + O (07 log? 6) , 
1— M% 
(7 = 6; K = a) 
K uy 2 Oe Ie ~y=(yt+1uny, ’ (5 2a) 
bn ty =O, (5. 2b) 
K ty + me as (vy + 1) (uy Ue), » (5. 3a) 
Ug = Vag = 05 (5. 3b) 


al: 
kK Us. = 5 Ug. se (y + 1) (uy %3) » 


F (5. 4a) 
2y—1) (y+ 
+ 04 (M+) —2y Ky my + co ia ) Ut U4; 
ee = Os (5. 4b) 
SS = 1+ 0(6® log?6) . (5.5) 


That is, the flow is irrotational in the first three orders. Potentials Pi Pa Ps 
can be introduced: 


1 . 
K ice ay Pipe + va Visa (y = 1) 7. ola’ ) (5.6) 


nee 
K G2, + Pape t+ Pap = Y +1) (G1, Y2,)a> (5.7) 


1 
K P34, + Pape F 3a = (vy + 1) (Pi, Ps,)x LY Piz Pry 
(5.8) 


ZW Ih +1 
—2yK m1, %,,+ (2 be ) 


2 
i Naga 


The boundary conditions on the body for these y’s are replaced by the 
expansion near the axis as carried out in equation (4.29). Further terms in 
the asymptotic expansions can be found as outlined in the previous section 
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and the results are 
a(%, #5 K) = Sy(x) log? + g(x; K) + Flog??)(2E* si sy) ] 
¥ — S t N" r Ol 
ely Seer Ly + FtSr ei) = 2 (y+ 1) S$, SY — K St] 
y2 73 eer oe 
+ [FOF (Si SY) = & +2) (Sel) + KS ne) 
+ (y +1) (Get) — gf K] + 0G log?) 


as r->Onearthe body 0<x<1, 


Ly : : Dy 6. ae 3 “ 
Po(x, 7; K) = S_(x) logy + g(x; K) + te (S; S3)' 7? log? + O(r? logy) , 


(5.10) 
93(*, 7; K) = S, S; log?7 + O(log?) (5.11) 
where the source strengths S,, S, are known from the body shape 
3 (a) = F(x) F(x), (5.12a) 
So(x) = —2 S, S) = —2 FF'(FF’)’ (5.12b) 


and the g,, g. are unknown functions. In general, to find g,, g, the differential 
equations must be solved and as will be seen in more detail later these functions 
are important for the pressure distribution. It is to be noted that in the expan- 
sion on the boundary the largest terms due to , and q, are of the same order. 
However it can be shown that the largest term due to “4, v4 in the expansion on 
the boundary is of higher order than those due to @, gs so that the expansion 
can be broken off. 

The differential equation (5.6) for g, is the classical transonic equation. It 
is non-linear and of changing type. The equation for , is however linear with 
variable coefficients depending on g,. The equation for ys and all the succeeding 
equations are also linear with variable coefficients and they may have various 
known forcing terms on the right hand side. 


6. Surface Pressures, Drag and Similarity 
The calculation of surface pressures depends on the velocity components 
on the boundary 7 = 62F. Using the expansions (5.9), (5.10), and (5.11), 


together with the expressions for the source strengths (5.12) an expansion for 


ZAMP VIII/2 
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the velocity components on the boundary can be written 


Gal% OF) _ 14 §2logd (2 Si) + 62 (Si logF + gi) + O(64log?6),  (6.1a) 


Wi 
Esa — 6F’ + 0(63 logs) . (6.1b) 
Next, the relation 

BS Ty (ol aoe 6.2 
ee Te ( : ) + 0(6?) (6.2) 

which holds throughout is solved for P/P,,: 
eee Rs 9 \We- 2) Py 6.3 
py = (gr lew (aR0(64 (6.3) 


An expansion of a?/U? on the boundary can be found from the invariant (3.4) 
by using (6.1): 


2 yd 

yx = 1— 6? logd[2 (y — 1) Si] 
(6.4) 
Zt: 


+ *(K —(y—1) [S;logF + 1 + sr]) + 0(64 log?6) . 


Thus, on the boundary 


P i | é Si | 4 2 
Py = 1 dogs (2y Si) — d?y(S; logF + gi + pe) + O(d4log?d) (6.5) 
and 
Beet 

oe — ae = — 6? logé (4 S”) 

Uy ge FMB, _ (6.6) 

4 ut 3 Monts 
— 622 (s; logF + gi + ir) + 0(54 log?6) . 


Equation (6.6) expresses the pressure on the body in terms of the body 
shape function F(x), a source strength S, = FF’, and an unknown function 
g(x). g1(x) has to be found from the solution of a boundary value problem for 
gi. Thus g, depends on K, g, = g,(x; K). The dominant term O(6? log 6) in 
(6.6) is part of the contribution of the sources along the axis (the other part is - 
6? Sj logF). A corresponding term due to sources exists in linearized supersonic 
theory where the expansion starts ([4], p. 89) 


Cp = — 6? logo (2 S}) + 0(62) (linearized subsonic or supersonic). (6.7) 


Thus, the sources contribute a dominant term in transonic flow which is twice 
as large as that in supersonic or subsonic flow. From a practical point of view, 
however, the term O(6?) is appreciable in comparison with the first term. 
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The term g,(x; K) represents the effect of the non-linear interaction be- 
tween sources and is the essential transonic part of the problem. The result is 
emphasized by a first-order similarity rule for the pressure coefficient. Define 
a quantity Cp is 


~ S2 
Cp =r + 2 Sj log (6? F) + 21 = —2 gi(x; K) + 0(6? log?). (6.8) 


The first order similarity rule is 


G. = Cp (x;K) with an error O(6?log?6). (6.9) 


The rule states that for bodies with similar shapes, the effect of the sources 
can be subtracted out so that there is a transonic similarity for the non-linear 
interaction effects. The rule in a form equivalent to (6.8) was given by Oswa- 
TITSCH and BERNDT [5]. An equivalent form of the rule can be written as 
Cp. > Se + 4 Sj) logé = f,(x; K) + O(6? log?6) . (6.10) 
Corresponding to the similarity law for Cp given by equation (6.9) a simi- 
larity law for the drag coefficient can be given. Defining 


2 767 
Cp = aor dx | (6.11) 


where 


; 2 : : zo 
A,, = dimensionless maximum cross-sectional area = 


it is found that for a body with a pointed nose and blunt base 


€, = SP + 8 FAL) F’%(1) log[o? F(1)] = —16 [ FF’ ¢; dx = j,(K) + 0(6% log?) 
0 


(6.12a) 
or 


af 
Cy=-8 {Cp FF dex. (6.12b) 
é 


The pressure on the base has arbitrarily been put equal to P,, in this derivation. 


7. Drag of Bodies of Nearly Circular Cross-Section 


The calculation of drag may be carried out in an analogous manner for a 
body whose cross-section is not quite circular. The surface of such a body will 
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be represented by 
S(«, 1, 0) =0=1r— 6F(x) V1 + tG(x, 8) (7.1) 


where t <1 and 6 is again the thickness ratio of the body of revolution 
vy = OF (x). If G(x, 0) is chosen so that 


2 G(x, 0) dd =0 (7.2) 


then the distribution of cross-sectional area is the same as for the body of 


revolution. In particular, the maximum cross-section area Ay = 2 67/4, and — 


this relation might be used as the definition of 6. The parameter 7 is considered 
to be independent of 6, and no assumption is made concerning the relative 
magnitudes of the two quantities. 

Since it has been shown that the flow over a body of revolution is approxi- 
mately irrotational [equation (5.5)], it would be expected that the initial 
terms of the velocity expansions for the present case may also be represented 
by means of potential functions. A new radial coordinate r* = 7/6 is introduced, 
and the following form of expansion is assumed: 


D(x, 7, 0; M.,,, 0, 7) 
UL 


= x+ 6? 9, (x, 7; K) + 64logé , (x, 7; K) +--- 
> Pylt, 0) or (4,7*, 07K) (7.3) 
+ fa (T, 0) pt (x, 7*, 0; K) ++. 


where the functions ,(x,7; K) are the solutions, assumed known, for t = 0. 
The boundary condition at the surface is obtained from equation (3.13a) if 


F(x, 0) is replaced by F(x) Vi+t G(x, 0). Substituting the expansion (7.3) 


d 


O° H1,(x, eFV1+7G) + & logd Pox (%, 6° FV1+4 7G) SS 


ret 
5 [FF m Tht Gal 


FYl+tG 


[1+ 0 oy. (2, PF VI+7G) +--+ uy ot (x FVI LG, 0) +++] 


G ene ADCS 
Oe, Ia a + TGR [v4 V7, (x, FV1 PitG, 0) + ve]. 


+ 4 @t, (% FVI+7G, 0) + 42 of (x,FV147G, 6) +--- : 


| 
t 
| 
| 
| 
| 
; 
| 
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Replacing the functions ,(x, 7) by their known expansions near the axis, and 
performing Taylor expansions about r* = F 


, 


ey 1 3 ( 
83} (1-5 1G+50Gt+ +) 4 5®logd “9291 (1-5 G+ +) fives 


FG 
+ a [ett F, 0) + =) px. .(x, F, 6) + | + #2 [ox (a, F, 0) +++] 


Lyx pe 


=; [FFP +57F@),] (1-F1G+2cG+...) 


x (1+ 62 logd 2S, +4+--:) 


23 Go 


ee se (+ +++) [oa 9%, ap a as cata 


Terms independent of t cancel because of the definition of y;(x, 7). The largest 
remaining terms determine mw, and the boundary condition for g¥*(x, 7*, 6): 


=O", (7. 6a) 
ot, (%, F, 9) = a (F*G),. (7.6b) 


The expansion has been chosen so that boundary conditions are applied 
on the mean surface F. Now, depending on the magnitude of T, “, will be either 
t 64log6 or t? 62. In the discussion of equation (7.16) it will be pointed out that 
terms linear in t cannot contribute to the drag, so we will take 


[lg = T? 0, (7. 7a) 


G FG G 
gt .(x, F, 8) = — = (F*G),— >" of. (4 FO) + 5% ot (x, F, 6) (7.70) 


without necessarily requiring that O(r 64 logd) < O(t? 6”). 

The differential equation to be satisfied is essentially the same as equa- 
tion (4.4): 
[1+ Ket ---—(y-N OM, +:-—Y-1 tH oe +: 

i 
ES [K 3 (oe 6* (M1,; = - %;) = eg 2 Kr od? iL, 
1 1 
ant (ot... + se ¥i,. + ar vt,,) + | 


: (7.8) 
= [+e mt +18 ot t+) ae 


1 0 1 0 
BAO Gye Ae OCT.) 6 * OF + pes (of, +--+) 4] 


x (Og bee $7 GT ++) | 
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The equations for p*(x, 7*, 0) are obtained by requiring that (7.8) be satisfied 
to every order of magnitude in the neighborhood of the body. It can be seen 
that the choice of 4, and is not influenced by the differential equation, and 
that w* and gy must satisfy 


u 1 F 
Cee UL. ar ge UE ee (7.9) 


fash Hed 190 


Expanding G(x, 0) in a Fourier series, 


= 3 G(x) e”® (7.10) 


n=—C 


where G(x) = 0 and G_, (x) = complex conjugate of G,(x). Substitution of 
this expression in the boundary condition (7.6b) allows solution for p¥(x,7*, 6). 


ot, (4 F, 0) = sie" BAG, of, | 
rare a (7.11) 
1 me rin 
py (*, 7*, 8) 2 Je, || vin Os : | 
The boundary condition (7.7b) becomes 
gx, (%, F, 6) = als a G, xe D>, |nl & Gy) e*°) 
co co n es a 
Fl. OES e! eth dal (F2G,)'e " 
and the solution for m#(x, 7*, 0) is 
eae ee i Bia. 
oh (a7, 0)=—7 DY Cla) a x Aa ein (7.12) 
where C;,(x) is defined by 
C, (x) eth 
k=—oo 


-» Se (<G, | 20 | GV eee ae fer Galea) 


-0o mM=—0o 


with C(x) taken equal to zero so that the solution dies out as 7* + co. 
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The pressure is obtained by expanding equation (6.3), with the use of the 
invariant (3.4). At the body surface the result is 


Cp = —2[ 68 gy (x, OF V1 + 7G)4+~— + 18% ot (x, FV1+76, 0) | 
+126? ot (x, FV14+7G, 6) +--+] 


[0% g(x, 8&FVI4+7G) $-. 4782 gt, (%, FV1+7G, 6) } (7.13) 


+ 1? 6? pF, (x, FV1+ =. 0) +...) 


; ve 
PP acy oe (He FV + 2G, 8) + +P +e 


i, 1S obtained by differentiating (5.9), and the terms representing the radial 


velocity component are replaced by the right-hand side of the boundary 
condition, in the form given by equation (7.5). After all terms have been 
expanded in Taylor series about 7* = F, 


Cp = -2 (89S; flog F)+ 5 7G-FeGt+-| ar ete, 


| 
+10 |g (x F, + Zt Got («FF 8) +++] 
: +2262 [pF (x, F, 6) ++] 4-)) 


—(2 [Fes frre) [1-5 res eats] (7.14) 


x [1 + 67logd6 2S, +---] 


Gite : : 
tox ght ee ota F 6) ++) 
ee AC ee : 


A considerable simplification occurs in the expression for drag. The drag 
coefficient is defined by a generalization of equation (6.11): 


MA 


1 27 . ce 
Does [dx [Cp oF V1 + 2G 4. (OF V1 + 2G) a8, 
pare (7.15) 


276 ‘ 
Sok fae [FF + 7 (F*6),| 40. 
A Gerth g 
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Substituting for Cp and combining terms when possible, 


1 2% 
Coa Gea t [ dx [ Cp, (F*G), 20 
0 0 


i 27 
= a | dx | [F* (F2G), + FF" F*G + 2FF' 9? (x, F, 0) d0 
0 0 


it 22 ¥ 
4 i 1 2 oS 
+ = 125? / dx / _1 FF (3 F° + FF") G 
Ss a | : L (7.16) 


— | F*(4F" + FF") GG 


x 


_ 3 FF G2— (F*G), oh (2, F, 8) 


— F* FG o* (x, F, 0) —2F F' gf (x, F, 6) 
2 . eae. do 
— F* G, of (x, F, 0) — 5 (a, F, 6)| 


where Cp, and Cp, are the results obtained for the body of revolution r = dF (x). 

It can now be argued that all terms linear in t, including those neglected 
above which are of order t 64 logdé and smaller, will drop out upon integration 
over 6. A more elaborate discussion is possible, but it seems sufficient to note 
that throughout the preceding development the Fourier expansions of terms 
linear in t never contain a term which is independent of 6. It follows that there 
is no term in Cp which is linear in t, and no generality has been lost by taking 
Mz = T? 6? in equation (7. 7a). 

The last integral can be simplified by using relations of the form 


[ord=2a Sc, ae 


nm=—COo 


Each integral of a double sum may therefore be replaced by a single summation ’ 
to be integrated over x. After rearranging terms, it is found that the remaining 
integration can be carried out explicitly to give 


Cp = Gy +2 7262 Y 
; Sie (7.17) 
[Tar (PAG) (PSG. ,) BA G3 FG ae G!,)} tt ole | 


where it has been assumed that the body has a pointed nose. It can be seen 


Vol. VIII, 1957 Expansion Procedures and Similarity Laws Za) 


that the summation reduces to zero in three eases. For a body with a pointed 
base, F (1) = 0 and consequently all terms in the sum disappear. If instead the 
base is blunt, but the rate of change of shape with x is zero, then 


Fp =G(L 0) = 0, 
and the sum again vanishes. Lastly, the sum can be zero if 
G(1, 6) = G,(t, 0) = 0, 


i. e., if the base is blunt and circular and if the derivative of the shape pertur- 
bation G(x, 9) is zero at the base. 

For the three types of bodies described, therefore, the change in drag due 
to small deviations from a circular cross-section is of smaller order than 1? 62: 


Cy = Cp, + 0(t? 62) (7.18) 


where 6 represents the order of magnitude of the body thickness ratio, and t 
the order of the perturbations in cross-section shape. This statement is actually 
a type of ‘area rule’, since the non-circular body is compared to a body of 
revolution having everywhere the same cross-sectional area. 

It is felt however that the analysis is not sufficiently general to allow 
application of the results to wing-body combinations. A separate expansion 
procedure is probably needed to give a good representation of a wing. 
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Zusammenfassung 


Die schallnahe Naherungslésung fiir die Umstrémung eines schlanken Dreh- 
kérpers wird mit Hilfe einer mathematischen Entwicklung beschrieben. Die 
Entwicklung wird durch eine abnehmende Folge von Funktionen des Dicken- 
verhaltnisses dargestellt, und ein Ahnlichkeitsparameter, der Machsche Zahl und 
Dickenverhaltnis enthalt, wird eingefiihrt. Die ersten Glieder der Entwicklung 
der Geschwindigkeiten nahe der Achse werden erhalten, und die Ahnlichkeits- 
sitze fiir Druckkoeffizient und Widerstandsbeiwert eines nichtangestellten Kor- 
pers werden abgeleitet. Mit Hilfe einer anderen Entwicklung wird eine angend- 
herte «area rule» (Flachenregel) fiir einen Kérper, dessen Querschnittsflache 
nicht kreisférmig ist, erhalten. 


(Received: June 25, 1956.) 
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Herstellung von Zufallszahlen auf Rechenautomaten 


Zusammenfassender Bericht 


Von SEBASTIAN VON HoERNER, GOttingen?) 


Zusammenfassung 


Zur Durchfithrung von Monte-Carlo-Methoden auf schnellen Maschinen ist 
die laufende Erzeugung langer Folgen von Zufallszahlen wahrend der Rech- 
nung notig. Es kommen rechnerische und physikalische Herstellung in Frage. 


I. Rechnerische Herstellung. Einige vorhandene Methoden der Herstellung 
werden diskutiert. Die Begriffe: Durchmischung, feste Zyklen, stabile Zyklen 
und Degeneration werden definiert und untersucht. 

Es wird eine neue Methode (Quadrieren und Erganzen) vorgeschlagen, die 
stabile Zyklen und Degeneration vermeidet und bei der die (stets unvermeid- 
lichen) festen Zyklen der Theorie folgen. Sie kann atch bei Maschinen geringer 
Stellenzahl allen praktischen Anforderungen geniigen. Ein ausfiithrlicher Test 
soll in Kitirze folgen. 


II. Physikalische Herstellung. Als Beispiel wird die Abzahlung von Gamma- 
Quanten diskutiert. Jeder Zahlrohrimpuls lasst einen Flipflop kippen, dessen 
Stellung als Zufallsduale tiber ein Gate abgefragt wird. — Zur Herstellung einer 
Folge von 108 Zufallsdualen geniigen im Mittel 5 Teilchen/Abfrage. 

Der Einfluss der Lange « der Abfrageimpulse und der Kippdauer y des 
Flipflop wird abgeschatzt. Um eine brauchbare Folge von rund 10® Zufalls- 
dualen mit einer Geschwindigkeit von 800 Dualen/s erzeugen zu kénnen, muss 
Oe esis Sell. 

Eine einfache Versuchsanordnung bestiitigte die Erwartungen und kénnte 
_ den Anforderungen der G2 geniigen. 


_Einleitung 


1. Unter der Bezeichnung « Monte-Carlo-Methode» werden seit einigen 
Jahren eine Reihe verschiedener Verfahren zusammengefasst, die in vielen 
Fallen eine ndherungsweise Behandlung sonst schwer zugdnglicher Differen- 
tialgleichungen und Integro-Differentialgleichungen erlauben: Das mathema- 
tische Problem wird auf einen Zufallsprozess abgebildet und dieser mit Hilfe 


1) Max-Planck-Institut fiir Physik. 
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von Zufallszahlen (in Art eines Wiirfelspieles) durchgefiihrt. Zur Einfiithrung 
empfehlen wir die zusammenfassenden Arbeiten [1] und paaay 

Die benétigten Zufallszahlen kann man durch Wiirfeln erhalten oder vor- 
handenen Tabellen («random numbers») entnehmen [3]. — Diese Moglichkeiten 
entfallen jedoch, wenn auf schnellen Maschinen gerechnet werden soll. Die Ein- 
gabe vieler Zufallszahlen wiirde meist weit linger dauern als die’ eigentliche 
Rechnung. Es ist also nétig, die Zahlen durch die Maschine selbst wahrend der 
Rechnung herstellen zu lassen. — Ganz allgemein bieten sich hierzu zwei Mog- 
lichkeiten an: rechnerische und physikalische Herstellung. 


2. Zur rechnerischen Herstellung miissen wir zunachst den Begriff der « Zu- 
fallszahl» naher betrachten. Denn jede Rechenvorschrift liefert festgelegte, 
determinierte Zahlenfolgen und nicht zufallige. Fiir die praktische Benutzung 
jedoch ist es ganz gleichgiiltig, ob diese Zahlen auf determinierte Weise ent- 
standen sind oder nicht. Wichtig ist nur, dass sie statistisch die gleichen Ergeb- 
nisse liefern wie ideale Zufallszahlen. Und da diese Verwendbarkeit sich testen 
lasst, so wollen wir einfach sagen: eine Folge von Zahlen nennen wir dann 
brauchbare Zufallszahlen, wenn sie eine Reihe bestimmter Tests erfiillt. Das 
bedeutet so etwa, dass die Folge keine groben Gesetzmassigkeiten erkennen 
lassen darf. Die Determiniertheit ihrer Herstellung darf man ihr nicht «an- 
sehen» kénnen. 

Wir werden zeigen, dass man in der Praxis keine beliebig langen Folgen 
herstellen kann. Die sich notwendig einstellenden Zyklen kénnen jedoch so 
gross gemacht werden, dass sie die Anwendung noch nicht stéren. 


3. Als phystkalische Herstellung werden wir in Abschnitt II die Abzahlung 
physikalischer Elementarprozesse benutzen. Wegen deren statistischer Natur 
sieht dies zundchst unproblematischer aus; es zeigt sich aber, dass man auch 
hierbei keine idealen Zufallsfolgen erhalten kann. Man kann jedoch Bedin- 
gungen angeben, unter denen die Folgen genitigend gut sind. 

Der Vorteil der physikalischen Herstellung liegt erstens in der grdsseren 
Geschwindigkeit : die Herstellung einer Zufallszahl ist nur ein Maschinenbefehl, 
wahrend fiir die rechnerische Herstellung je nach Verfahren 10—20 Befehle 
nétig sind. Und zweitens ist die Problematik der so erzeugten Zahlen zweifellos 
geringer und besser zu tiberblicken als bei rechnerischer Herstellung. Dagegen 
hat die rechnerische Herstellung den Vorteil, keine zusatzliche Apparatur zu 
benétigen. Und zweitens liefert sie reproduzierbare Folgen, was fiir die Fehler- 
suche und die Méglichkeit der Wiederholung ganz wesentlich ist (vgl. NEU- 
MANN [1)). ' 


4, Weiterhin ist zu tiberlegen, ob man getestete Folgen benutzen will oder 
nicht. Wir wollen hierunter Folgen verstehen, die auf Grund bestimmter Tests 


2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 52. 
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aus einer grésseren Anzahl von Folgen ausgewahlt worden sind, wahrend der 
Rest verworfen wurde. 

Bei ungetesteten Folgen ist das Ergebnis einer Rechnung in vollem Umfang 
dem Zufall ausgesetzt. Man erhalt eine «natiirliche Streuung», die im allge- 
meinen ebenso gross sein wird wie bei Ausfiihrung des entsprechenden Expe- 
rimentes (zum Beispiel Messung einer Diffusionskonstanten durch Abzahlung 
von Teilchen). ; 

Die Benutzung getesteter Folgen dagegen kann die Streuung je nach Zu- 
sammenpassen von Tests und Anwendung unter Umstanden stark reduzieren. 
Andererseits ist zu bedenken, dass der Erwartungswert E einer Funktion f 
einer Zufallsvariablen x nicht gleich der Funktion des Erwartungswertes von 
x ist, sondern angendhert : 

d*f 


E(H)} = HEC) + 5} (FE) neo" 


Kommen in der Rechnung nichtlineare Funktionen der Zufallszahlen vor, so 
wird eine reduzierte Streuung o im allgemeinen das Ergebnis verfalschen. 


5. Schliesslich ist noch zu fragen, wie der Test einer Methode durchzufihren 
ist. Die meisten in der Literatur vorhandenen Beispiele gehen dabei so vor: die 
zu testende Methode wird so angewandt, dass sie bei einer langeren Folge von 
Zahlen noch keine Fehler zeigen sollte. Und dies wird durch eine Reihe spe- 
zieller Tests ( (Haufigkeiten von Ziffern, Ziffernkombinationen, Liicken zwi- 
schen Ziffern usw.) gepriift. Mit Hilfe der 7?-Methode entscheidet man dann, 
ob die Hypothese « Zufallszahl » anzunehmen oder zu verwerfen ist. 

Eine andere Méglichkeit ist folgende. Die Methode wird so angewandt, dass 
sie bei einer nicht zu langen Folge bereits ihre Fehler zeigt. Ausserdem wird 
die Art und Grésse der Fehler theoretisch abgeschatzt. Der Vergleich zwischen 
Theorie und Beobachtung zeigt dann, ob die Fehler der Methode durch die 
Theorie richtig erfasst worden sind. Ist das der Fall, so lassen sich die Bedin- 
gungen angeben, unter denen die Methode brauchbare Ergebnisse liefern wird. 

Zweifellos wird die erste Art der Priifung stets angewandt werden miissen. — 
Wir méchten aber hiermit einigen Nachdruck auf die zweite Art der Priifung 
legen. Wir meinen, dass die direkte Berechnung und Messung von Fehlern 
mehr Vertrauen gibt als die blosse Feststellung, dass die Fehler bei einer - 
Beobachtung « unter der Messgenauigkeit » gelegen haben. 


I. Rechnerische Herstellung 
A. Vorbemerkung 


1. Zunachst kénnte man daran denken, zum Beispiel die Folge der Dezimal- 
ziffern einer transzendenten Zahl als Zufallsfolge zu benutzen. Bei dieser und 
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allen ahnlichen Methoden ist jedoch klar, dass mit verniinftigem Rechenauf- 
wand sich keine langeren Folgen herstellen lassen. 

Deshalb wird man praktisch wohl ganz allgemein so vorgehen miissen, dass 
man durch Anwendung einer speziellen (kurzen) Rechenoperation R auf die 
letzten mw Zufallszahlen die nachste Zufallszahl erhalt. Zur Verminderung des 
Aufwandes sollte w klein sein, als einfachsten Fall wird man, falls modglich, 
f= 1 wahlen: 


X41 = R(x). (1) 


2. Dabei wird es (je nach der Wahl von R) oft so sein, dass eine bestimmte 
Ziffer (Nummer &) von x;,,; nur von einem Teil (Anzahl d) der Ziffern von x; 
abhangt, 

Meare = Rs Kieegt ree s Berga)» (2) 
wahrend die tibrigen Ziffern von x; zur betrachteten Ziffer von x;,, keinen 
Beitrag liefern. Wir wollen d den Grad der Durchmischung nennen, der im all- 
gemeinen von & abhangen wird: d(k), und der héchstens gleich der Anzahl s der 
benutzten Ziffern sein kann: d Ss... 

Je kleiner die Durchmischung ist, um so eher wird man der Folge x; grobe 
Gesetzmassigkeiten ansehen k6énnen. Ein Beispiel liefert Abschnitt C. 


3. Wir wollen nun der Reihe nach einige Méglichkeiten fiir die Wahl von R 
durchgehen. Eine der einfachsten Operationen ware die Addition einer festen 
oder laufenden Zahl. Dabei ist aber die Durchmischung sehr klein, zum Beispiel 
fiir die letzte Ziffer, k = 1, ist d(1) =1, und im Mittel ist (wegen Ubertragung 
aus den riickwartigen Ziffern) etwa d = 1,5. Die Addition entfallt also wegen 
zu geringer Durchmischung. 


4. Als nachste Méglichkeit kénnte man an die Multiplikation mit einer 
festen Zahl denken, unter Weiterbenutzung nur der letzten s Ziffern. Dies 
fiihrt zu den Methoden von Topp, MosHMAN und Juncosa von Abschnitt C. 

Beginnen wir die Zahlung von k mit der letzten Stelle, so ist hierbei 
d(k) =k. Die Durchmischung ist also gut bis auf eine Anzahl der letzten 
Ziffern. 


5. Bei der Multiplikation zweier gleichlanger Zahlen kénnen wir eine noch 
héhere Durchmischung erreichen, wenn wir nicht die letzten s Ziffern benutzen, 
sondern die mittleren s Ziffern. Dann ist 


: 2S 
5 ixss und a=. (3) 


Das einfachste Beispiel hierfiir ist das « Quadrieren und Abschneiden» nach 
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von, NEumANN (Abschnitt D). Dies hat jedoch den Nachteil, dass die Folgen 
«degenerieren» kénnen, das heisst in eine Folge von Nullen ausarten. 


6. Eine Methode, welche dies vermeidet, wird in Abschnitt E vorgeschlagen 
und diskutiert. 


B. Allgemeine Fehlerquellen 
a) Feste Zyklen 


1. Da jede Maschine nur mit einer endlichen Stellenzahl s rechnet, muss 
jede Rechenvorschrift der Art von (1) notwendig zu Zyklen fiihren — ganz 
gleich, welche spezielle Form man fir (1) wahlt. Denn bei s Ziffern ist der 
Wertevorrat N (Anzahl der ttberhaupt méglichen Zahlen) : 


N= 10''dezimal bzw. ‘N= 2° dual. (4) 


Nach spatestens N Zahlen muss ich somit eine Zahl erhalten, die schon einmal 
da war. Dann miissen aber wegen der Eindeutigkeit von (1) sich auch die 
folgenden Zahlen alle gleichen; ich erhalte somit einen Zyklus der Lange 
P<=N,. Diese (allein durch die endliche Stellenzahl s der Maschine bedingten) 
Zyklen wollen wir feste Zyklen nennen. 


2. Im folgenden setzen wir eine Rechenvorschrift (1) voraus, die ausser 
diesen unvermeidlichen festen Zyklen sonst keinerlei Gesetzmassigkeiten er- 
kennen lasst. (Zum Beispiel kénnte fiir jede der N Zahlen x; ein Nachfolger 
%;,.1 durch Wirfeln festgelegt worden sein.) Dies wollen wir eine ideale Methode 
nennen. 

Als erstes stellen wir die Frage: wie gross ist (bei einer idealen Methode) 
die Wahrscheinlichkeit /(), dass man (von einer beliebigen Zahl ausgehend) 
genau  Zahlen erzeugt, bevor man erstmalig zu einer Zahl kommt, die schon 
einmal da war? Denn nur diese Folge von  Zahlen ware als Zufallsfolge 
brauchbar. — Diese Wahrscheinlichkeit f(m) ist das Produkt der Wahrschein- 
lichkeiten dafiir, dass ich die ersten Zahlen ohne Wiederholung erhalte, und 
dafiir, dass die Zahl Nummer 1 + 1 eine schon vorhandene ist: 


fee ina) (rep lce eee 


Indem wir (5) logarithmieren, unter der Voraussetzung N > n entwickeln 
und nach Potenzen von 1/N sortieren, erhalten wir 


Hn) = se "VEN (n8 <6 WN). (6) 
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Setzen wir noch ) N > 1 voraus, so liefert die Verteilung (6): 


Maximum = 75 = VN, 


Mittelwert = m = Ve N= 125 VN, 


It 


Streuung = o = (2 = ) NN = 0,655 VN, 


2 


——— 
——: 


mittlere Periodenlinge = P = ae | 
Die Formeln (6) und (7) kénnen wir benutzen, um eine gegebene Rechenvor- 
schrift (1) darauf zu priifen, wie gut sie in bezug auf die festen Zyklen einer 
idealen Methode entspricht. Man nimmt mit kleiner Stellenzahl s alle N még- 
lichen Zahlen als Ausgang einer Folge und notiert. jeweils die Anzahl 1 der 
ohne Wiederholung erzeugten Zahlen. Abschnitt E zeigt ein Beispiel. 


3. Weiterhin fragen wir nach der Wahrscheinlichkeit F(m) dafiir, dass wir 
uns bei einer Folge von » Zahlen bereits in einem Zyklus befinden: 


F(n) = / f(n) dn= —" (F€1). (8) 


Formel (8) benutzen wir, um zu entscheiden, wieviel Zahlen wir mit einer 
gegebenen Maschine (bei idealer Methode) erzeugen kénnen, ohne mit Zyklen 
rechnen zu miissen. Dabei muss fiir F eine geeignete Sicherheitsgrenze will- 
kiirlich festgelegt werden: wir halten F = 10-% fiir ausreichend (jede tausendste 
Rechnung kénnte Fehler durch Zyklen aufweisen). 

Hat zum Beispiel eine Maschine nur sechs Dezimalen, so erhalten wir mit 
F==-10-*% aus. (8) 8, 
n = 2000 ~ 45 
als Grenze fiir die Lange brauchbarer Folgen. Das heisst aber, dass Maschinen 
geringer Stellenzahl zur Herstellung von Zufallszahlen nicht geeignet sind. — 
Eine Méglichkeit der Verbesserung besteht in der Benutzung nicht nur der 
letzten Zahl zur Erzeugung der nichsten, also zum Beispiel = 2 oder 3 

(siehe Abschnitt E unter f). 

Als weiteres Beispiel diene die Géttinger Maschine G2. Mit rund 30 Be- 
fehlen/s, einem Programm von 10h und rund 20 Befehlen fiir Herstellung und 
Verarbeitung einer Zufallszahl wiirde eine Folge von = 5 - 10* Zufallszahlen 
_bendtigt. Mit s = 15 Dezimalen ist dann 


F=10-°, 


“was eine véllig ausreichende Sicherheit bedeutet. 
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b) Nicht herstellbare Zahlen 


Bei einer idealen Rechenvorschrift ist nicht jede Zahl aus einer anderen 
herstellbar. Veranschaulichen wir wieder die ideale Methode durch Auswirfeln 
eines Nachfolgers fiir jede der N Zahlen, so ist die Wahrscheinlichkeit dafir, 
dass eine herausgegriffene Zahl nicht als Nachfolger einer der N Zahlen ge- 
wiirfelt worden ist: Hor 

(1 — x) wen, (9) 


N 


Fiir die mittlere Anzahl w der nicht herstellbaren Zahlen ergab eine langere 


Rechnun a 
e u = 0,413 N + 0,338 VN. (10) 


Die knappe Halfte aller Zahlen sind also nicht herstellbar. Der ganzen Zahl : 


darf man bei der Verwendung somit nie eine spezielle Bedeutung zulegen, 
sondern stets nur einer kleineren Anzahl von Ziffern. Im tibrigen kann auch 
(10) bei kleiner Stellenzahl zum Testen einer Methode benutzt werden. 


c) Stabile Zyklen 


Auch bei beliebig grosser Stellenzahl der Maschine ist das Auftreten von 
Zyklen prinzipiell méglich. Zum Beispiel ergibt jede Zahl x, die der Gleichung 


Rigas (11) 
gentigt, einen Zyklus der Periodenlange P = 1. 


Wir stellen nun die Frage: Kann ein solcher Zyklus «sich einspielen »? Wir 
beginnen eine Folge in der Nahe eines Zyklus 


Rw + &) = x+e (&<1). 


Ist nun ¢, kleiner als €9, so wollen wir die Lésung von (11) einen stabilen Zyklus — 


nennen. In gleicher Weise definiert 


EES ee: (12) | 


einen Zyklus der Lange P, der dann stabil ist, wenn aus der Form von Ry 
folgt, dass 


4 


Ejipo Ss 6; (13) | 


ist. Im allgemeinen wird dies nur dann der Fall sein, wenn wir uns schon 
gentigend nahe an dem Zyklus befinden: 


€; = Ein ° (14) 


Diese obere Grenze ¢,, wollen wir den Stabilitdtsbereich des Zyklus nennen. 
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Bei einer idealen Methode sollte es keine stabilen Zyklen geben. Brauchbar 


ist eine Methode nur dann, wenn die Sumime ihrer Stabilitatsbereiche hin- 
reichend klein ist. 


d) Degeneration 


Unter Degeneration wollen wir, in Ubereinstimmung mit [1], den Fall ver- 
stehen, dass in bestimmten Positionen einer Zahlenfolge sich die Ziffern repro- 
duzieren, wahrend der Rest weiter variiert. Ein bekanntes Beispiel liefert die 
v. Neumannsche Methode, wenn die zweite Halfte einer Zahl aus Nullen 
besteht. Der Grund hierfiir ist, dass die Null bei der Multiplikation eine aus- 
gezeichnete Rolle spielt. Das Produkt ist unabhangig vom zweiten Faktor, die 
Durchmischung ist reduziert. 


C. Multiplikation mit fester Zahl 


1. Die Herstellung der nachsten Zufallszahl aus der vorherigen durch Mul- 
tiplikation mit einer festen Zahl A, unter Benutzung der jeweils letzten s 
Ziffern des Produktes, schreiben wir 


%.4, = R(x) =A-x, (mod 10°). (15) 


Der von der Position k abhangige Grad der Durchmischung wurde in A 4. mit 
d(k) =k angegeben. Dies hat zur Folge, dass die Periodenlange der festen 
Zyklen von k abhangt. 

Betrachten wir zunachst die letzte Ziffer (k = 1) von x;,,. Sie steht allein 
unter dem Einfluss der letzten Ziffer von x; und der (festen) letzten Ziffer von A. 
Wir kénnen sie somit unabhangig von den anderen Ziffern betrachten. [hr 
Wertevorrat ist N = 10, wir miissen also feste Zyklen erhalten mit P, < 10. 
In gleicher Weise kann man dann zeigen, dass die Ziffern der Position k 
periodisch sein miissen mit 

et ra =. 10% (16) 


Dies gilt fiir beliebige, aber feste Zahlen A. Stabile Zyklen kénnen mit A el 
nicht auftreten, da ¢;,, = A ¢; ist. 


2. Will man die Lange der Zyklen direkt wissen an Stelle der Unglei- 
chung (16), so muss A in spezieller Weise gewahlt werden. Dies fihrt zu den 
Methoden [4] von J. MosHman und M.L. Juncosa fir dezimale Maschinen 
und von O. Taussky fiir duale Maschinen. Ausfiihrliche Tests ergaben eine gute 
Brauchbarkeit der vorderen Ziffern, wahrend die (sechs) letzten Dezimalen 
_ wegen zu geringer Periodenlange entfielen. Betreffs Einzelheiten verweisen wir 
auf die Originalarbeit. 


_ ZAMP VIII/3 
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3. Vielleicht sollten zwei Dinge noch untersucht werden. 
a) Angenommen, die Maschine macht einen Rechenfehler. Was geschieht 


von da ab? 

b) Es sollten auch Ziffernkombinationen getestet werden (« serial tests»). 
In [1] berichtet ForsyTue, dass eine der vy. Neumannschen verwandte Methode 
wegen dieses Testes verworfen werden musste, obwohl alle anderen Tests gute 
Resultate gaben. 

Im iibrigen soll nochmals darauf hingewiesen werden, dass die von MosH- 
MAN angegebenen Periodenlangen sich nur auf die vorderste Ziffer beziehen. 
Fiir die folgenden Ziffern ist in seinen Formeln der Wert von s um je 1 pro 
Stelle zu erniedrigen. 


D. Quadrieren und Abschneiden (v. NEUMANN{1]) 


a) Methode 


Die letzte s-ziffrige Zahl wird quadriert, und die mittleren s Ziffern des 
Produktes werden als neue Zahl benutzt. Ist zum Beispiel s = 4 und die letzte 
Zahl %;= 2573, $0 ist 47 = 6617750 und 47.4 = Ola 


b) Stabile Zyklen und Degeneration 


1. Wir nehmen an, es gabe einen Zyklus der Lange P nach (12). Stabil 
kann er nur sein, wenn (13) erfiillt ist. Durch P-malige Anwendung unserer 
Operation R lautet (13) bis auf Glieder héherer Ordnung jetzt: 


(2a D092) (2,4 LO) Sul Oia ce aa ee (17) 


Dies Produkt der P Faktoren kann aber nur dann < 1 sein, wenn mindestens — 
ein Faktor < 1 ist. Er habe den Index 7. Es muss also gelten 


x; 05 0 ee 
Ist dies der Fall, dann ist aber laut Herstellungsvorschrift 
Bic 020° 10°? usw. 


Das heisst, die Folge lauft nach Null, ohne sich vorher wiederholen zu kénnen. 
Die Voraussetzung des Zyklus war also falsch, falls nicht alle x =O sind. 
Somit ist « = O der einzige stabile Zyklus der Methode. 

Der Stabilitatsbereich dieses Zyklus ergibt sich aus 


Rie) == 6.0 n> ere 
zu 
e, = 10%, (18) 
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Besteht die vordere Hilfte einer Zahl aus Nullen, so lauft die Folge von da ab 
in den Null-Zyklus ein. 


2. Besteht die hintere Halfte einer Zahl aus Nullen, so hat jede folgende 
Zahl die gleiche Eigenschaft, wahrend die vordere Halfte weiter variiert. Die 
Folge ist degeneriert. 


3. Die Wahrscheinlichkeit dafiir, dass eine herausgegriffene Zahl der Folge 


in der ersten oder der zweiten Hialfte (oder tiberhaupt) nur Nullen hat, betragt 
fiir grosse NV 


agree 

VN 

Damit betragt die Wahrscheinlichkeit Aj(m) dafiir, dass man wahrend einer 
Folge der Lange bereits in eine solche Nullenfolge eingelaufen ist 


Z2ai\h 2H 
his 1S (1 cm. : ain (<1). (19) 
Vergleichen wir dies mit der entsprechenden Wahrscheinlichkeit F fiir feste 
Zyklen (8) F = n?/2 N, so ist 

72 
Fe}. (20) 
8 

und das bedeutet, dass bei dieser Methode die festen Zyklen gegeniiber den 
Nullenfolgen keine Rolle spielen. 


4. Schliesslich wurde noch die Wahrscheinlichkeit g(m) berechnet, genau 
beim (mz + 1)-ten Schritt eine solche Nullenfolge zu beginnen, 


8") = (1 ae) ap (21) 
und ihr Mittelwert m, 
my = 3" n a(n) =5VN (m1). (22) 


5. Nehmen wir als Beispiel wieder die G2 mit N = 10”, so gibt (19) fiir 
eine Folge der Lange = 5- 10* eine Sicherheit von Fy = 0,003; die Wahr- 
scheinlichkeit, durch Nullenfolgen Fehler zu erhalten, ist also eigentlich 


bereits zu gross. 
Dagegen hatten wir nach (22) im Mittel erst nach m, = 1,5: 10’ Zahlen 


den Beginn einer Nullenfolge zu erwarten. 
c) Ziffernhadufigkert 


1. In [1] weist P. C. Hammer darauf hin, dass die so erzeugten Zahlen nicht 
genau einer Gleichverteilung entsprechen, denn falls x gleichverteilt ist, so 


36 SEBASTIAN VON HOERNER ZAMP 


gehorcht y = x? der Verteilung y~'?/2. Dies gibt eine Bevorzugung kleiner 
Zahlen, die sich in den herausgeschnittenen s mittleren Dezimalen allerdings 
nur noch so wenig auswirkt, dass eine Folge von tiber 3000 zehnziffrigen 
Zahlen in [1] eine Reihe verschiedener Tests gut bestanden hat. 


2. Wir wollen nun den hieraus entstehenden Fehler theoretisch abschatzen. 
Die Zahlen x und y denken wir uns als Dezimalbruch mit s bzw. 2s Ziffern 
geschrieben in den Grenzen 0... 1. Dann ist, wie gesagt, 

1 


iy) = Sie (23) 


Wir fragen nun nach der Verteilung g,(D), wobei D die k-te Dezimale hinter 
dem Komma sein soll. Zum Beispiel ist fiir die erste Dezimale 


1/10 
D WES ome pley 
g,(D) | tao + #) # as VD+1-YD), (24) 
und ahnlich giit allgemein 
jok-1_1 107% 
2(D)= > ? fv 10'-* + D-10-* 4.2) dz. (25) 
v=0 


Eine Bevorzugung kleiner Ziffern wird sich am starksten bei der Null zeigen. 
Wir fragen deshalb nur nach der Haufigkeit der Null in der k-ten Ziffer hinter 
dem Komma. Durch (23) ist f gegeben, wir integrieren und erhalten 


Lae 


g,(0) = 10-"2 (V10»+1—V10»). (26) 


0 


Zur Berechnung der Summe rechneten wir die ersten drei Glieder direkt aus 

und ersetzten den Rest durch ein Integral. Nach einiger Zwischenrechnung 

ergibt sich ee 
147,61 or 

i (27) 


Bei Gleichverteilung wirden wir 1/10 erhalten. Als relativen Fehler 4 wollen - 
wir bezeichnen (1/10) 
uch Sod Sie a8 - 10-42 
A, 1/10 = Ole ae (28) 
3. Um zu entscheiden, ob dieser Fehler zulassig ist oder nicht, vergleichen 
wir ihn mit der Streuung. Bei Gleichverteilung wiirden wir bei Herstellung von 
m Zahlen in jeder Position im Mittel m= 1/10 Nullen zu erwarten haben, mit 


einer Streuung von o = (3/10) Vn. Wir wollen den Fehler der Methode dann 
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wieder als zulassig betrachten, wenn er kleiner als 1 /3 der relativen Streuung ist: 


ae: 


1 
A er ee ee, 
3 m yn 


(29) 
Ist s wieder die Anzahl der benutzten Ziffern, so hat die erste Ziffer die Position 
k= 1+ s/2, und Forderung (29) lautet 


n<0,17VN (30) 


als héchste Lange einer brauchbaren Folge. 


d) Diskussion 


Betrachten wir zu einem Vergleich die Formeln (19), (20) und (30), so sehen 
wir, dass die praktische Begrenzung der Methode allein durch die Wahrschein- 
lichkeit -) der Nullenfolgen gegeben ist. Fordern wir wieder F< 10-* als 
ausreichende Sicherheit, so ist die maximale Lange einer brauchbaren Folge 
von s-ziffrigen Zahlen 


n<5-10-*VN =5- 108-4, (31) 


Das bedeutet aber, dass zur Erzeugung langerer Folgen nur Maschinen sehr 
hoher Stellenzahl in Frage kommen. Auch bei 14 Dezimalziffern dirften nur 
Folgen von héchstens n = 5000 Zahlen benutzt werden. 

In Abschnitt Ba3. wiesen wir darauf hin, dass die Wahrscheinlichkeit fester 
Zyklen sich verringern lasst, wenn zur Erzeugung der nachsten Zahl mehrere 
der letzten Zahlen benutzt werden. Das gleiche gilt auch fiir die Nullenfolgen. 
Zwei Beispiele einer Anderung der Methode gibt ForsyTHE [1] an, zum Beispiel 
statt des Quadrates der letzten Zahl wird das Produkt der beiden letzten 
Zahlen genommen. 

Wir mochten aber lieber so vorgehen, dass wir zunachst eine Methode suchen, 
die auch bei Benutzung nur der letzten Zahl den Anforderungen schon geniigt. 
Die Benutzung mehrerer Zahlen wird dann erst bei Maschinen extrem geringer 
Stellenzahl nétig. 


E. Quadrieren und Erganzen 


a) Methode 


Wir suchen also eine Methode, die 
1. nur die letzte Zahl benutzt, das heisst Form (1) hat, 
. méglichst einfach ist, 
. keine stabilen Zyklen hat, 
. keine Degeneration aufweist, 
. relativ gut zu iiberblicken ist, 


ak wWhd 
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und méchten dafiir vorschlagen: 


# eee (G2 Oe | 
¥41 = R(x,) = | (32) 
\ta (<2). | 
Alle x sollen in 6 < «<1 liegen. — Von dieser Folge wird die erste Ziffer 


sicher keine brauchbare Zufallsreihe sein, die beiden nachsten wohl auch noch 
nicht. Wir werden aber abschitzen, von welcher Ziffer ab der Einfluss der 
Herstellung zu vernachlassigen ist. 


b) Allgemeine Bedingungen 


Damit das Verfahren verniinftig lauft, miissen wir zwei Bedingungen fiir 
die Konstanten stellen. Die erzeugten Zahlen sollen erstens unterhalb von 1 
bleiben: 

a+o<1, (33) 
und zweitens oberhalb von b: . 


be+a> ob. (34) 


In Figur 1 darf der Arbeitspunkt (a, 6) also nur zwischen der Geraden (33) und 
der Parabel (34) liegen, am besten méglichst weit von beiden Linien entfernt. 
Die so erzeugten Zahlen liegen dann in 


bSx<b+a. (35) 


Degeneration oder ahnliches sind bei dieser Methode nicht méglich. 


c) Stabile Zyklen 


Wir setzen nach (12) einen Zyklus der Lange P voraus und beginnen in 
seiner Nahe eine Zahlenfolge nach (32); nach P-maliger Anwendung von (32) 
lautet Bedingung (13) dann bis auf Glieder hdherer Ordnung 


ERO ER oT “ise ee ee | (36) 


L 


Der Zyklus kann nur stabil sein, wenn (36) erfiillt ist. Dies lasst sich nun ganz 
leicht vermeiden, indem wir durch 


uh 
D7 ae (37) 


dafiir sorgen, dass nach (35) alle x grésser als 1/2 sind. Dann ist kein stabiler 
Zyklus moglich. 
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In Figur 1 muss also unser Arbeitspunkt (a, 6) oberhalb der gestrichelien 
Linie liegen, am besten in der Nahe des eingetragenen Punktes. (Wiirden wir 
dagegen zum Beispiel a = 0,6 und 6 = 0,1 wahlen, so wiirde jede Folge nach 
kurzer Zeit in einen stabilen Zyklus der Lange 3 einlaufen.) 


10 


Figur 1 


Wahl des Arbeitspunktes (a, b) fiir die Herstellungsvorschrift (82). Oberhalb von 6 = 1/2 ist kein 
stabiler Zyklus méglich. Bester Arbeitspunkt etwa bei @. 


Figur 2 
jufigkei S h genau ” Schritten. 
Haufigkeit des Auftretens fester Zyklen nac Z 
— theoretische Verteilung nach (6); s Stellenzahl (dezimal); N= a- LO* = Wer tevorrat; gs 
rischer Mittelwert von (7); 1 theoretischer Mittelwert nach (7); P empirische mittlere Perioden- 
lange (theor. = mo/2). 
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d) Feste Zyklen 


Wir wollen priifen, ob diese Methode beziiglich ihrer festen Zyklen sich wie 
eine «ideale Methode» von Abschnitt Ba) verhalt. Nennen wir » wieder die 
Linge einer Folge bis zur ersten Wiederholung einer bereits vorhandenen 
Zahl, so kénnen wir fiir geringe Stellenzahl s jede der N méglichen Zahlen als 
Anfangszahl einer Folge wahlen und  bestimmen. Da alle Zahlen zwischen 
b...6+ a liegen, ist N =a- 10°. 

Dieser Test der Methode wurde durchgefiihrt mit s = 2 und s = 3. Die er- 
haltenen Verteilungen von n zeigt Figur 2, die durchgezogenen Linien sind 
die theoretischen Verteilungen nach (6). Wir sehen, dass sowohl m als auch P 
wie auch die Form der Verteilung in beiden Fallen sehr schon mit der Theorie 
iibereinstimmen. Beziiglich ihrer festen Zyklen scheint also unsere Methode den 
Anforderungen der idealen Methode zu entsprechen. 


e) Ziffernhaufigkert 


1. Ahnlich wie bei der v. Neumannschen Methode sind auch hier die er- 
zeugten Zahlen nicht gleich verteilt. Um abzuschatzen, von welcher Dezimale 
(oder Duale) an wir den Fehler vernachlassigen kénnen, miissten wir die Ver- 
teilung f(x) kennen, die sich bei langen Folgen einspielt. Analytisch ist sie 
dadurch definiert, dass sie durch Anwendung der Operation (32) auf das 
Argument ihre Gestalt nicht andert. 

Die Funktion /(x) ist nicht differenzierbar. Da /(x) an den Randern x = db 
und *=6+ a Spriinge (der ungefahren Grdsse 1/a) hat, und da f(x) sich 
durch Anwendung von (32) auf x nicht andern soll, so muss /(x) auch an den 
Stellen R(b) = b?+ a und R(b+ a) = (b+ a)? Spriinge besitzen. Und von 
diesen beiden Unstetigkeiten kénnen wir auf gleiche Weise auf zwei weitere 
Unstetigkeiten schliessen usw. Dabei vermindert sich die Sprungstarke jeweils 
um einen Faktor 7 = dx/dR = 1/2 x, wobei x hier der Ort der vorigen Unstetig- 
keit ist. Im Mittel ist (in Umgebung des vorgeschlagenen Arbeitspunktes) 


1 
Wegen (35) und (37) ist stets 7 <1, wir erhalten somit eine abzahlbar - 
unendliche Folge von Unstetigkeiten abnehmender Starke. Dabei gehorcht 
die Dichte der Unstetigkeiten der gleichen Verteilung { wie die Zahlen x. Fiir 


die folgende Abschatzung kéunen wir jedoch die Verteilung der Unstetigkeiten 
Isa Gleichverteilung ansetzen. 


2. Wir fragen nun nach der Abweichung der Verteilung der k-ten Dermat 
von einer Gleichverteilung. Da f(x) nicht differenzierbar ist, berechnen wir als 
erstes die quadratisch gemittelte Differenz # der Funktion f(x) an zwei benach- 
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barten Orten x und x + ¢, als Funktion von ce. —Im Intervall ¢ liegt ein Sprung 
der ersten Sorte (Starke nia) mit der Wahrscheinlichkeit 2 ¢/a. Mit gleicher 
Wahrscheinlichkeit und davon unabhingig liegt im Intervall noch ein Sprung 


der zweiten Sorte (Starke 7?/a), desgleichen der dritten Sorte, und so fort. 
Daraus folgt 


w= (fa) — e+ P= "| (2) + (LY + (ZV 4+-- 


a 


Normieren wir noch durch die mittlere Hohe 1/a der Verteilung f(x), so erhalten 
wir fiir die relative Differenz d(e) benachbarter /-Werte: 

h Ze 
(Ware / differenzierbar, so miisste d~e sein.) Um (40) zu priifen, wurde d 
fiir einige « empirisch bestimmt aus einer nach (32) hergestellten Folge von 
1300 Zahlen’). Die Ubereinstimmung mit (40) ist gut: 


| € | d nach (40) | d empirisch 
0,01 | 0,254 Q729 
0,02 0,360 0,40 
0,05 


0,567 | 0,56 


3. Die mittlere Differenz zweier Ziffern ist 5. Mit e =5-10-* und a= 0,3 


erhalten wir 
P57 lots (41) 


Innerhalb des Bereiches } ... b + a kommt nun eine bestimmte Ziffer der k-ten 
Dezimale a-10*~! mal vor. Die relative Haufigkeitsdifferenz A zweier Ziffern 
ist also die Summe iiber die a-10*~! einzelnen d-Werte, die bei statistischer 
Verteilung der Unstetigkeiten voneinander unabhangig sind. Wir erhalten 


somit im Mittel i 


—k 
TaD aie 10cke (42) 


Nennen wir g,(D) die Verteilung der Ziffern D in der k-ten Dezimale hinter dem 
Komma und ist |D, — D,| = 5, so ist A dabei wie folgt definiert : 


V(gn(Ds) — ex (D2)? 
ye “TO : (43) 


3) Fiir die Durchfithrung danke ich Herrn FissER. 
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4, Zur Kontrolle zeigt Figur 3 die Verteilung der zweiten Dezimale bei einer 
Folge von 1270 Zahlen. Daraus erhalten wir empirisch A = 0,41, was mit dem 
aus (42) folgenden Wert 0,33 im Einklang ist. Ftir die dritte Dezimale sollte aus 
(42) A = 0,033 sein. Figur 3 zeigt, dass der Effekt hier bereits in der Streuung 
untergeht, die bei Gleichverteilung 3/Vn = 0,084 (das Zweieinhalbfache) 
betriige. 


LA. 

Or oe ee beGe mors 

j= D 
Figur 3 


0123456789 


Die Haufigkeitsverteilung der 2. und 3. Dezimale bei einer Folge von n = 1270 Zahlen. 
Die gestrichelten Geraden sind Mittelwert - Streuung bei Gleichverteilung. 


Nun fordern wir wieder, dass dieser Methodenfehler héchstens 1/3 der 
Streuung betragen soll. Fiir die maximale Lange einer Folge, beziiglich ihrer 
k-ten Dezimale hinter dem Komma, erhalten wir damit: 


n < 102-3 (44) 


5. Nun war aber die Lange einer Folge bereits durch die Wahrscheinlichkeit 
F begrenzt (8), in einen festen Zyklus einzulaufen. Wir wollen deshalb die 
Frage stellen: von der wievielten Dezimale ab sind die Ziffern einer Folge 
brauchbar, deren Lange durch (8) gegeben ist ? (44) ist also jetzt eine Gleichung 
fiir k, wobei m durch (8) gegeben ist. 

Setzen wir wieder F = 10-3, a = 0,3, N = a- 10°, so erhalten wir schliesslich 
ungefahr 


beg (45) 


als Position der ersten brauchbaren Ziffer. Das heisst: stellen wir nach der 


Methode (32) eine Zahlenfolge maximaler Lange (8) her, so ist von jeder dieser 7 


Zahlen das erste Viertel der Ziffern zu streichen, der Rest ist als Zufallszahl 
benutzbar. 


f) Mittelung 


1. Normalerweise wird die durch (8) dargestellte Grenze m durchaus geniigen. 
Anders jedoch bei einer Maschine geringer Stellenzahl. Die Grenze lage so 
etwa bei 10 Dezimalstellen; dann ist » = 4000, und da wir nach (45) sieben 
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Ziffern benutzen diirfen, so gibt das rund 3 104 Ziffern insgesamt. Und dies 
diirfte ftir die meisten Probleme fiir ein « Spiel geniigen. 

Bei einer Maschine kleinerer Stellenzahl gehen wir so vor, dass wir statt der 
Form (1) die w letzten Zahlen zur-Herstellung der nichsten Zahl benutzen. — 
Der einfachste Fall ist die Mitteilung. Zunachst zeigen wir, dass wir yw auch 
wieder nicht zu gross wahlen diirfen. Denn der Mittelwert m; der letzten yu 
Zahlen unterscheidet sich fiir grosse w nur wenig von m;,,,, und folglich wiirden 
auch die hieraus erzeugten Zahlen x;,, und x;,5 dicht beieinander liegen. — 
Wir meinen aber, dass dieser Effekt zum Beispiel bei fu = 3 oder 4 noch keine 
Rolle spielen diirfte, da wir ja ohnehin die ersten Ziffern nicht als Zufallszahlen 
benutzen. Unsere Methode lautet also jetzt, mit zum Beispiel u = 3: 


ym; (m; 2 b), 
4,4 = Rin y= Ri (46) 
\m? +a (m2? < 5b) 
mit 
1 
Mm, = = (4%, %_1 + %_9) - 


Die drei Bedingungen (33), (34) und (37) beziiglich des Arbeitspunktes (a, 0) 
bleiben unverandert, desgleichen der Zusammenhang (44) zwischen der Lange 
einer Folge und der Anzahl benutzbarer Dezimalen. 


2. Neu zu berechnen ist die Wahrscheinlichkeit fester Zyklen. Um in einen 
festen Zyklus einzulaufen, ist jetzt ndtig, dass die letzten drei Zahlen gleich 
einem schon vorhandenen Zahlentripel sind. Die Anzahl vorhandener Tripel 
ist etwa gleich der Anzahl vorhandener Zahlen; und die Anzahl iiberhaupt 
méglicher Tripel (Wertevorrat) ist gleich der dritten Potenz des Wertevorrates 
a-10° der einzelnen Zahlen. Beides zusammen heisst, dass wir die Formeln 
des Abschnittes Ba) iitbernehmen kénnen, wenn wir 


Wee 10 a3 10" (47) 


einsetzen. Fordern wir wieder F = 10-%, so ergibt sich die maximale Lange 


einer Folge jetzt aus (8) zu 
os (8) nm =0,8 510°=* (48) 


Und fiir eine Folge dieser Lange gibt (44) dann 
k= ist+1 (49) 


als Position der ersten brauchbaren Dezimale. Wegen der grésseren Lange der 
Folgen ist also nur noch das letzte Viertel der Ziffern als Zufallszahl zu be- 
nutzen. Aus (49) séhen wir auch, dass es keinen Sinn hatte, w grésser als 4 zu 
wahlen. Denn von da an ist die maximale Lange einer Folge unabhangig von mu 
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und allein durch (44) mit k = s gegeben, da & nicht grosser als s gewahlt 
werden kann. 


3. Die endgiiltige Grenze der Methode liegt also bei u = 4 und alleiniger 
Benutzung der letzten Dezimalen als Zufallszahl. Die maximale Lange ist 


dann nach (44) Pee (50) 


was selbst bei extrem geringer Stellenzahl noch ausreichen diirfte (zum Beispiel 
selbst s = 4 gibt noch m = 10°). 


g) Diskussion 


Die Ergebnisse der letzten Abschnitte zeigen, dass die durch (32) vorge- 
schlagene Methode alle zu Beginn gestellten Forderungen gut erfillt. Bei 
richtig gewahltem Arbeitspunkt sind weder stabile Zyklen noch Degeneration 
méglich. Und die stets vorhandenen festen Zyklen entsprechen in ihrer Zahl 
und Lange denen einer «idealen » Methode. 

Der Einfluss der Nicht-Gleichverteilung wurde berechnet: mit Ausnahme 
des ersten Viertels der Ziffern ist er zu vernachlassigen. 

Fiir Maschinen mit mehr als 10 Dezimalen geniigt die Herstellung der nach- 
sten Zufallszahl aus der letzten allein. Fiir Maschinen mit geringerer Stellenzahl 
wird die Methode durch Mittelung iiber die letzten drei bis vier Zahlen nach 
(46) erweitert und geniigt dann allen praktisch vorkommenden Anspriichen. 

Ein ausfihrliches Testprogramm fiir langere Zahlenfolgen ist fiir die G2 
in Vorbereitung. Eine Arbeit von FISSER wird dariiber berichten. 


II. Physikalische Herstellung 
A. Vorbemerkung 


1. Auf eine allgemeine Behandlung der physikalischen Herstellung wollen 
wir hier verzichten. Wir méchten eine, anscheinend brauchbare, Methode vor- 
schlagen, deren Eigenschaften wir in Abschnitt B theoretisch abschatzen. In 
Abschnitt C beschreiben wir dann eine einfache Versuchsanordnung und die 
mit ihr erzielten Ergebnisse. 


2. Eine handliche Méglichkeit zur schnellen physikalischen Erzeugung von 
Zufallsdualen besteht in der Abzaéhlung von Teilchen eines radioaktiven Zer- 
falles mit Hilfe eines Zahlrohres (zur Durchfithrung wurden die Gammaquanten 
eines Radiumprdparates benutzt). Figur 4 zeigt das Schema einer solchen 
Anordnung. Jeder Impuls des Zahlrohres lasst einen Flipflop kippen, dessen 
einer Anodenausgang zu dem einen Steuergitter eines Gate fiihrt. Auf das 
zweite Steuergitter kommen aus der Rechenmaschine die Abfrageimpulse, die 
das Gate nur dann passieren kénnen, wenn der Flipflop-Ausgang positiv ist. 
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Es wird also registriert, ob im Augenblick der Abfrage die bisherige Gesamt- 
zahl der Teilchen gerade oder ungerade ist. Und je nach dem wird eine Eins 
oder Null in den Akkumulator eingeschrieben. 


| 
Zusatz-Apparatur Masohine 

| 

! 


“ __ u . 
= Zahlrohr Flipflop 


| 
Abfrage-Impulse 


Figur 4 
Allgemeines Schaltschema fiir Abzahlung von Gammaquanten. 


3. Beziiglich der zu stellenden Anforderungen haben wir uns nach der Gét- 
tinger Maschine G2 gerichtet. Alle Abschatzungen sind jedoch ganz allgemein 
gehalten und kénnen auch anderweitig verwendet werden. 

Es soll ein Maschinenbefehl geniigen, um wahrend einer Trommeldrehung 
den Akkumulator mit einer geniigenden Anzahl von Zufallsdualen zu fiillen. — 
Eine Trommeldrehung der G2 dauert 20 ms, und wahrend dieser Zeit finden 
16 Umlaufe des Akkumulators statt. Bei jedem Umlauf soll der Inhalt des 
Akkumulators mit 2 multipliziert werden, und eine Zufallsduale soll nach- 
gefiillt werden. Es miissen also 16 Dualen in 20 ms erzeugt werden, oder 
800 Dualen/s. 


B. Theoretische Abschatzung 
a) Anzahl der Teilchen|Abfrage 


1. Waren zum Beispiel die Abfragen haufiger als die Teilchen, so wirden 
immer mehrere aufeinanderfolgende Abfragen das gleiche Ergebnis liefern ; 
wir bekimen somit eine starke positive Korrelation benachbarter Dualen. 
Lassen wir nun die Frequenz der Abfragen fest und erhéhen die Teilchenzahl, 
so wird diese Korrelation kleiner. Unsere Frage lautet: Von wann ab ist sie zu 
vernachlassigen ? — Zunachst einige Definitionen: 


yu = Anzahl der Abfragen/s, 

t = 1/u = Zeitintervall zwischen zwei Abfragen, 

y =mittlere Teilchenzahl/s, (51) 
y = Teilchenzahl in herausgegriffenem Intervall t , 

m = v/u = mittlere Teilchenzahl in 7, 

n = Lange einer Folge von Zufallsdualen. 


46 SEBASTIAN VON HOERNER ZAMP 


Die Wahrscheinlichkeit f(y) dafiir, dass in einem herausgegriffenen Intervall 
der Lange t sich y Teilchen befinden, gentigt einer Poisson-Verteilung 


iy) = re (52) 
Damit ist die Wahrscheinlichkeit G (bzw. U) dafiir, dass y gerade (ungerade) ist 


G=(1+ 4-4 Prat eo feo EU Silas 4 te) eo a (53) 


Wir definieren als relativen Methodenfehler der Korrelation 


(G- — (Gf 
A, = 24a =2(E=0) (54) 
und erhalten 


2. Die Forderung, dass die Korrelation benachbarter Zahlen zu vernach- 
lassigen sei, heisst nun: der Methodenfehler A, soll einen bestimmten Bruchteil 
der bei Gleichverteilung vorhandenen Streuung nicht tiberschreiten, wobei wir 
1/3 fiir ausreichend halten wollen. 

Bei statistischer Gleichverteilung von geraden und ungeraden y ware die 
Streuung der durch (54) definierten Grésse, bei einer Folge von » Abfragen, 


1/V/n. Unsere Forderung lautet somit 


1 


Aue ae (56) 


und unter Benutzung dekadischer Logarithmen mit (55) : 


n m m > 0,895 + 0,576 log E (57) 
pf 3,2 Die nebenstehende Tabelle zeigt einige Mindestwerte von m 
a Fe fiir vorgegebene 7. Wir sehen daraus, dass mit m= 5 Teil- 
BF re chen/Abfrage bereits allen praktischen Anforderungen ge- 


168 a6 ntigt ist. (Auch wenn m um 10° schwanken sollte, so ist 
n = 10° noch zulassig.) 


b) Endliche Impulslange 


1. Eine weitere Fehlerquelle besteht darin, dass sowohl die Lange der Ab- 
frageimpulse als auch die Kippzeit des Flipflop nicht beliebig klein gemacht 
werden kénnen. Wir miissen also untersuchen, welchen Einfluss auf die Zif- 
fernhaufigkeit diejenigen Teilchen ergeben, die wahrend einer Abfrage erschei- 
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nen. Dabei hangt es vom Kathodenpotential des Gate ab, welcher Bruchteil 
der « abgeschnittenen » Abfrageimpulse noch als Eins gezahlt wird. 

‘Das Kathodenpotential V kann innerhalb eines Bereiches von } Volt 
gewahlt werden, dessen untere (linke) Grenze 7 dadurch gegeben ist, dass bei 
noch kleinerem Potential alle Abfrageimpulse das Gate passieren kénnen, 
ganz gleich, ob der Flipflop auf Eins oder Null steht. Und jenseits der oberen 
(rechten) Grenze r= 1+ 6 ist das Potential so hoch, dass tiberhaupt keine 
Abfrageimpulse durch das Gate gelangen. 


% 


100 
p 
50 
0 | 
Figur 5 
Haufigkeit der Eins. 
V Kathodenpotential der Gate-RGhre; p Prozentsatz der erzeugten Einsen; — —~—— verschwindende 


Impulslange; endliche Impulslinge. 

Tragen wir die relative Anzahl # der Einsen in Prozent tiber V auf, so 
hatten wir im Idealfalle verschwindender Impulsliangen und Kippzeiten: 
eae 00%, fir V<_t, p= 50% fir 1 < V <7,wnd = 0% fir x < V., Dies 
entspricht der gestrichelten Linie in Figur 5. 

Statt dessen haben wir bei endlichen Zeiten einen Verlauf nach Art der 
ausgezogenen Linie zu erwarten. Ihre Form hangt ab von der Form der Impulse 
und des Kippens sowie von der Charakteristik der Gate-Réhre. Wir wollen 
dies nicht naher untersuchen, sondern den Verlauf #(V) in erster Naherung als 
Gerade betrachten mit den Grenzwerten #; und #,. Thre negative Steigung ist 
n = (pi — f,)/b. Vo ist der zu wahlende Arbeitspunkt der Kathodenspannung, 
der aber nur fiir #, < 50% existiert. 


2. In Figur 6 schematisieren wir die Impulsformen zur Abschatzung von 

p, und #,. — Befinden wir uns an der linken Grenze J, so geniigt bereits eine 
verschwindend kleine Erhéhung des Flipflop-Ausganges, um eine Eins zu er 
zeugen. Damit die Abfrage eine Null ergibt, muss erstens ¢g < ¢, sein, das heisst, 
za einer Zeit < t, —y muss ein Teilchen begonnen haben, den Flipflop zu 
schliessen. Und zweitens darf ein Offnen des Flipflop nur mit ¢; > ¢, geschehen, 
das heisst, erst zu einer Zeit > ¢, darf ein Teilchen beginnen, den Flipflop zu 
 6ffnen. Das bedeutet aber, dass jedes in das Intervall ¢,— (4,—y)=a+y 
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| | = 

| | | | oe 
| | 1\ | 4 Ni 
/ | \ / \ 
t; a SET ITE 
aan yt ty t, by nc Gade 7 ' 

i ang Are we, 
Figur 6 


Schematisierte Impulsformen. 
a) Abfrage-Impulse; b) Flipflop-Ausgang. (Die auslésenden Teilchen liegen bei tf; und fz.) 


fallende Teilchen eine Eins erzeugt, gleichgiiltig, ob das Teilchen den Flipflop 
gedffnet oder geschlossen hat. 

Die Anzahl der auf diese Weise erzeugten Einsen/s betragt v mw («+ y). 
Ohne diesen Effekt sollten die Halfte davon Nullen sein. Es werden somit 


i 
7 Ve (at+y) 
zusatzliche Einsen/s erzeugt. Driicken wir nun noch das Verhdltnis dieser 


Grésse gegeniiber der Anzahl w der Abfragen/s in Prozent aus, so erhalten wir 
die Abweichung gq, an der linken Grenze 


1, =, — 50% =» (un +7) 50%. (58) 


Durch eine 4hnliche Ableitung erhalten wir an der rechten Grenze 
Und schliesslich folgt hieraus fiir die (negative) Steigung 7 


wo OPS. Pra Pn eo (Oe Poe % 
ai eae 3 b 20 rates (60) 


3. Aus (59) ersehen wir, dass y > 8 sein muss. Nur wenn dies erfiillt ist, 
gibt es einen Arbeitspunkt Vj, an dem die Einsen ebenso haufig sind wie die 
Nullen. Die Abfrageimpulse miissen also méglichst sfztz sein, zumindest muss . 
ihr konstanter Teil kiirzer sein als die Kippzeit des Flipflop. 


4. Nun kann aber die Spannung V, nicht beliebig genau eingestellt und 
konstant gehalten werden. Die Abweichung 


v=V,—V (61) 


mag etwa -+ 1 Volt betragen, um die Apparatur nicht unnétig zu komplizieren. 
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Der relative Methodenfehler der Ziffernhautigkeit ist damit 
ga ks gy 
Ay — She. b- (a =e B oF 2¥) : ‘ (62) 
Wir fordern nun wieder, dass dieser Methodenfehler kleiner sein soll als 


ein Drittel der relativen Streuung bei Gleichverteilung, ahnlich wie in (56). Da 
B <y sein muss, wollen wir f gegeniiber « + 2 y vernachlassigen : 


Ay =; @+2y) << —— (63) 
3)/n 
und somit 
pips wes \2 
sgn operg rer gr yine (64) 


Oder umgekehrt erhalten wir mit: 6 = 50 Volt, v= -++1 Volt, » = 4000/s, m= 108: 


a+ 2y<4,2-10-°s. (65) 


C. Versuchsanordnung 
a) Aufbau 


1. Zur Priifung der Ergebnisse des vorigen Abschnittes wurde eine einfache 
Apparatur gebaut*). Gegeniiber der prinzipiellen Schaltung der Figur 4 sind 
einige Zusdtze nétig. Erstens liegt zwischen Zahlrohr und Flipflop ein Verstar- 
ker. Zweitens zeigte sich, dass die lange Abfallszeit (100 us) der Zahlrohrimpulse 
durch Differenzieren auf knapp 3 us verkiirzt werden konnte. Drittens verhin- 
dert ein Diskriminator (Normierung der Impulshéhen), dass der Flipflop, auf 
Grund innerer Unsymmetrie, durch zu kleine Impulse nur noch in einer Rich- 
tung gekippt werden kann, was eine ungleiche Haufigkeit von Eins und Null 
ergeben wiirde. 

Das Praparat war nicht zu stark: in 1 m Entfernung konnte man sich noch 
dauernd aufhalten. Das Zahlrohr wurde mit 1,3 kV betrieben, und der Verstar- 
ker hatte einen Faktor von 104. Die Schaltung des Flipflop wurde einer Arbeit 
von Piroty jr. [5] entnommen, wobei der Eingang tiber eine Duodiode an beide 
Gitter gelegt wurde. 


2. Es standen zwei Zahler zur Verfiigung mit jeweils hintereinandergeschal- 
teten zwei dekadischen Zahlréhren und mechanischem Zahlwerk. Die Zahler 
waren bis etwa 1500 Impulsen/s brauchbar, und bei Hintereinanderschaltung 


aller 4 Zahlréhren bis etwa 25000/s. 


4) Fiir die Benutzung der benétigten Apparaturen méchte ich Herrn Prof. Wirtz danken, 
sowie Herrn BERNHARD Meyer fiir viele praktische Ratschlage. 
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Die Kippzeit des Flipflop (einschliesslich etwas langsamerem Einschwingen) 
betrug y = 1,5 us. Eine Zahlung der Impulse vor und hinter dem Flipflop 
ergab, dass er auch bei 15000 Teilchen/s noch bei jedem Teilchen kippte. 

Zur Abfrage wurden die Impulse des Oszillographen benutzt, deren kleinste 
Lange « = 4 us betrug. Die beiden Zahler waren zur Messung gleichzeitig an 
Abfrageimpulse und Gate-Ausgang gelegt. 


b) Ergebnis 


1. Es wurden fiinf Messreihen durchgefiihrt, eine davon (Nr. 2) zeigt Figur 7. 
Das gesamte Ergebnis ist dargestellt in der folgenden Tabelle: 


Nr. A 


ik 50000 4 0 | 25 | 4700 | 800} 1,3} -- 0,35] 0,8} — 0,3 6 | 0,142 
Z 10000 | 12 6 |) 75 |-4700)) 800 | 3,2 | 42951) 3,6) 405827 sOmi2o 
3 23000 | 28 | 25 | 65 | 4700 | 800} 6,9) + 5,5 | 6,6} + 3,5 | 23 | 0,135 
4 
5 


IZO00F pelZ, 6 | 75 | 4700 | 210} 3,2) + 1,1 3,4 | — 0,7 26 | 0,160 
250000 | 12 6 | 75 | 1000 | 210} 0,7} +0,2 | 08] — 0,2 26 | 0,035 


ee ee 
berechnet gemessen 


(58) und (59) 


A Anzahl der Abfragen/Messpunkt; a, B, h nach Figur 6; » Teilchen/s; 
uw Abfragen/s; g, 6, 7 nach Figur 5. 


tl 


es ees 
60 70 80 90 100 
= l Volt ——— 


Figur 7 - 
Messreihe q(V) Nr. 2, entsprechend Figur 5. 
Mita = 12s,“ = 800 Abfragen/s, » = 4700 Teilchen/s. Jeder Messpunkt stellt 10000 Abfragen dar. 
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2. Die Ubereinstimmung der berechneten und der gemessenen g,-Werte 
ist recht gut. Dagegen zeigen die gemessenen g,-Werte zwar allgemein, den 
gleichen Gang wie die berechneten, liegen aber zumeist niedriger. Der Grund 
dafiir ist, dass die Abfrageimpulse nicht die idealisierte Form der Figur 6 be- 
sitzen, sondern auch in ihrem waagerechten Teil 6 etwas gekritmmt sind, was 
durch Verzerrung innerhalb der Apparatur noch verstarkt wurde. Somit ist die 
konstante Lange f effektiv kleiner als in der Tabelle angegeben. 

Eine weitere Bestatigung unserer Abschatzung zeigt sich darin, dass die 
Steigung 7 durch Herabsetzen der Abfragefrequenz ju (Messreihe Nr. 4) nicht 
vermindert wurde, sondern erst durch Herabsetzen der Teilchenfrequenz » 
(Messreihe Nr. 5). 


c) Verwendbarkeit 


1. Nachdem unsere Abschatzungen durch die Versuchsapparatur bestatigt 
wurden, fragen wir nach dem Nutzeffekt, der sich mit unveranderter Apparatur 
bei einem Anschluss an die G2 ergeben wiirde. 

Die Anzahl der Abfragen muss, wie gesagt, “ = 800/s betragen, und die 
Teilchenzahl muss somit etwa v = 4000/s sein. — Die Abfrageimpulse kénnen 
der G2 entnommen werden mit einer Lange « = 2 us (f ~ 0) und einer Hohe 
h = 140 V, wodurch sich der Bereich } auf etwa 50 V vergréssern diirfte. 

Damit ist dann die maximale Lange einer Folge von Zufallsdualen nach 
(64) : 

n= 7 = 10°; (66) 


was fiir alle praktischen Zwecke stets ausreichen diirfte. Das heisst also, dass 
unsere Versuchsapparatur ohne prinzipielle Anderungen den Anforderungen 
der G2 bereits geniigen wiirde. 


2. Zum Abschluss fassen wir nochmals einige Bedingungen zusammen, die 

bei Verwendung der Methode an anderen Maschinen erfiillt sein miissen. Wir 

gehen davon aus, dass eine spezielle Maschine w Abfragen/s benotigt. Dagegen 

wollen wir 6 = 50 V und v= + 1 V als feste Zahlen in die Rechnung setzen. 

Die Apparatur soll brauchbare Folgen bis zu einer Lange von n= 10° Zufalls- 
dualen erzeugen. Mit diesen Werten erhalten wir aus (57) 


y=Su (67) 


fiir die mittlere Teilchenzahl/s. Und beziiglich der Lange « der Abfrageimpulse 
und der Kippdauer y des Flipflop erhalten wir mit (64) 


1 
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Die konstante Lange f der Abfrageimpulse muss kleiner sein als die Kippzeit 


des Flipflop 
a de (69) 


Sind diese drei Bedingungen erfiillt, so betragen die Methodenfehler héch- 
stens ein Drittel der Streuung, liegen also noch « unterhalb der Messgenauigkeit ». 
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Summary 


I. Some of the available methods for the arithmetical generation of random 
numbers on electronical computers are discussed. The basic notions (mixing of 
digits, fixed and stable cycles and degeneration) and different ways of testing are 
introduced and investigated. 

A new method is developed [see formula (32)] avoiding stable cycles and de- 
generation. Theoretical estimates and several tests gave full agreement. 


Il. As an example of the physical generation of random numbers, the counting 
of Gammaquanta is treated theoretically. The estimates have been confirmed 
using a simple physical generator, built in the computer G2 of the Max- Planck- 
Institut fiir Physik, Gdttingen. 


(Eingegangen: 28. Februar 1956.) 
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Circular Cylinder Under Periodic Fluctuations of Temperature 


By VActav VopicéKa, Plzen, Czechoslovakia!) 


Summary 


This paper brings general results about heat waves in homogeneous isotropic 
cylindrical bodies in periodic temperature fields. The following deductions are 
of technical significance and imply as simple special cases the results presented 
usually in literature*). With the purpose of giving to engineers and physicists 
reliable formulae, we take pains of detailed calculations. 


§ 1. Statement of the Problem 


The boundaries of the cylinder 0 < @ <7, |z| Ss are subjected to periodic 
changes of temperature, namely: the curved surface 0 = r to the fluctuations 
f(y, z) e*', the upper base z= s to g(e, gy) e*”’ and z= —s to ho, p) e*™"; 
f(y, 2), g(e, @), Ale, y) are given functions, wp , 1, @, constants, ¢ denotes the 


time and i = /—1. On separate parts of the surface we take into account the 
convection of heat, respectively with the coefficients Ay, h,, i, and have to 
determine the resulting steady temperature distribution u = u/(o, ¢, z; ¢) in the 
solid. 

The case in question is expressed mathematically (with usual notation) by 
the following equations: 


Ou 2 [| Ou 1 Ou if 0*u 07u 
AR aa ( 2 : ates Ree 2 ) , | 
ot 00" o 00 0 Op 0z 


C2657 Slap ate, \see sc, te 0. 


OH + glu — fp.) =0, 9-7. @) 
oH su elo, 9) 4 =0, 2=5: 
Zz 

@) 
24 _igiu—ho.g) ei] =0, z= 5. 


The whole problem can evidently be reduced to three simple cases by putting 
U = Uy + Uy + Uo; by, = 4,(0,9, 232), &=1, 2, 3. (4) 


1) Technische Hochschule. 
2) V. VopiéKa, Schweiz. Arch. angew. Wiss. 6 (1948); Appl. sci. Res. [A] 5 (1955). 
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Each of the functions , “1, % has to satisfy the equation (1), but the other 
conditions are 


1 Aube s Ou ao Oy 5: 
“ao + ho Lto — Fp, z)er t= 0, 0 =F; i +hyu=0, 2=S; 


(5) 
OU, thie eeaet 
wy, hints == OS s 
88 igus =0, o=7) “thin slog) e@]=0, z=5; 
do 0z (6) 
Ou 
— = hy Uy — 0) , 2=>-—S. 
ie het, =O, C= 7; eae 0, Sis 
(7) 
0 tw. 
2 — Ig [up — ho, y) e**4] =0, z=—s. 


The well-known Bernoulli method gives the following normal form of 
particular integrals v, = v,(0, y, 2;¢) of (1), suitable for solving our three 
partial problems: 


VA0, 0,2; 1) = Ro) DP io) Z(z) Tt) sta A ee (8.1) 
Herein 
R,(o) = 1 (= 0), @,(y) = A® cosng + A® sinny , 
n n\ » n n n (8 2) 
Oz 0z Aa*t. : 


Z(z) =cos—+Csin—, Tih= 4"; c= Vast+ 8; | 


6 real, A are free parameters. A), A®) and C constants of integration. 


§ 2. Determination of Uy 


We assume wu, to be linearly composed from our normal integrals (8.1). 
Substituting v, in the two last conditions (5) gives 


C (6d cosé + A, s sind) — d sind + h, scosd =0, 


C (d cosé + hg s sind) + 6 sind — hys cosd = 0, 
and from this we obtain on the one hand the equation 
(hy hy s* — 67) sin2 6 + (hy + hy) s dcos25=0, (10) 


on the other 
0 sind — hy s cosé 
‘beos dd by Ssind (11) 
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Therefore, we easily get by (8.2) 


R,,,(@) = 1,(=# 0) , D,,,(p) = AS COS @ + Ay i sin Pp ; | 
3 12 
Op C08 Oy (1 - =) + hys Sin 6, (1= a) (ee | me 
Anal == ~ 0, Cos 6, + hy s sind, ee Ne VA s* + 05 
and (8.1) yields a doubly infinite sequence 
0. G2; t) = KR, le). D,.(y) 2,02) oh"; on = 0/1, 2)00., R12 ae © 13) 


of particular solutions with the period 2 in @ and satisfying for all », k the 
two last conditions (5). The values 6, entering in (12) and (13) are positive 
roots of (10). 

Thus, all te calculations reduce to determining the parameter / and 
the coefficients A\'), A‘) in the expression 


aay 
Uo (0, YP, 2; t) = 3S yale 0,9, 2, t) a de ds Bunl@) y) 2,(2) (14) 
n=0 k= 


according to the first condition (5), 1. e. when primes mean derivatives 


det 57 3 Real) ) + lo Ryx(7)] Ona(P) 242) = ho Hp, 2) e** - 


n=0 k= 
This gives firstl , 
3 b Bae (15) 
a 
and secondly 
DS AR Raw”) + Fo Rua) Zp(2) = Ao fnal2)i = 9,1,2,..., (16-1) 
1 
> AP Rig”) + fo Ruel) Ze(2) = Mo fnal2)i m= 1,2,--. (16.2) 
k=1 
where 
1 ia eae COB 1 oe bee 
fosl@) = af Hes 2) do , fps fies BE gorse hy 2 OD) 
0 


But we have, according to the two last relations (5), for each pair of values 
Re 1 (k,l =1, 2, ...) swecessively 


s Ss 


ae Ui 
[40 EE Sars ee |e Zi — Z, 2) dz 
- s? y / 5 pe 
= gg nti %Alt--2= 9. 
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Consequently, any function g(z) satisfying the well-known conditions can be 
expressed in the form 


Ja (6) d 
g(z) = D7 G™ Z,(z), G® 7 ‘ (17) 
k=1 epee dt 
Thus, the equations (16.1)—(16.3) lead to 
h (k) 
oe Rox”) a hy Ro, (7) : a Rin ) ho Rap”) i (18.1) 
B= D2 ie kee pore ie 
where 
La “) Z,(0) at i fla, ¢) de , 
Fy EY Dee s ; 
is Z2(0) de 
/ nee -(18.2) 
Fe) , eae er m : 
Ee . iS 
cs [BUS a | 
Sarl 2, athe, Oo ae ky Cy ee 


Using (12), (14), (15), (18.1) and (18.2) gives finally the required solution 
Uy = U0, Y, 2; t) of our first partial problem in the form 


m= ae tS ae ae 5a [fet f(a, 


[50 dex 


Rox (@) (g) 
Ray +h, Re = a er ec wa co" a). 


, 
a 


The expressions R,,,(0) and Z;(z) are defined in (12). 
It is easy to accomplish the indicated calculations. Before all, we have 
from (12) and (15) at once 


Rial?) + hy R,, (7) = [(Ao r+ n) I (e,) + & I ERAC 


|/ez- : pee 


and so it remains to evaluate the integral in the denominator of (19). 


(19.1) 


il 
Yr 
Y 
we 
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On the one hand we find by means of (5) 


; Of 2 a) J ” 7. 
/ (+ 22 — 7,3) dz =| (—2, 2, = 4,7) dz=Z, Z, | = hy Z2(s) + hy Z2(—s) , 


on the other we obtain from (8.2) and (11) for 6 = 6, 


Adding these two results together gives 


S 


n ‘ 2 
| 23(2) dz = s (1+ CR + So thy ZH(s) + hy Z2(—5)] , 
“ ae 


= 


or again, by (11) and (12) 


S{2 Of (Of + hy s*) + s [hy df + hy (d, cos2 6, + hy s sin2 6,)*}} 


2 dj (0, cos d, + hy s sin d,)? (2) 


| HO a = 
§ 3. Finding of u, 


We again assume u, in form of a linear combination of particular integrals 
v,(0, p, z; #), as given by (8.1) and (8.2). 
Substituting v, in the first and the last equations (6) gives 


éJ,,,(€) + (or +) Le) =90, C(dcosd +h, ssind) + dsind—h,scosd=0 (20) 
and we obtain a doubly infinite set of particular integrals 
etre 2: 1) ==R(0) D. fe) Z da) eos tee O, 1,2, Rd, 2 ZL) 
with 
Ry (0) = T,{“** Q) , Pyal—) = Auk cosnp + Arh sin ny , 
(2122) 


by : & 
On 7, COS On (1+ <) + hy ssindy, (1+ =} . | 


’ 


Zq4(2) = 3,» C08 5, , + hy S SIND yy 
&,, are the roots of the first equation (20), 6,, and A are connected by the 
relation (8.2). The integrals (21.1) satisfy not only the basic equation (1), but 
also the first and the last relation (6). 
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Putting 


0, 9,431) = 3S male 9, 1:1) = 08 SY Rilo) (9) Zane) (22) 


n=0 k= n=O0k=1 


we find from the second condition (6) 
Jet SS Zi al6) + I Zyl) Ruel) Paul) = He (0, @) 2 


This leads to 


ite eee (23) 
a 
and 
DS An [Zn(s) + Ay Z,,(8)] Ryx(0) = 1 G10); = 0,1,2,..., (24.1) 
eal 
ora Zins) + Ay Zy2(8)] R,n(0) = M1 Snol); 2 =1,2,... (24.2) 
pal 
with 
27 
ee Oe &,1(0) cosua ,. 
Enle) = 397 | sle-a) dx , a) -+ fe oe ) iano" ; | (24.3) 
Hie b2, | 


Now, we have for each pair of indices k + 1 (k, 1 = 1, 2,...) by the theory 
of cylindrical functions (primes denote the o-derivatives) 


and from this in the well-known way ~ see also the first equation (20) — 


re fe Runle) Rule) do ={ 0[Ria(0) Rule) — Ryalo) Rewlell Yeo = 0- 
0 


Hence, we have obtained the following relations of orthogonality 


(25) 
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and each function g(o) of the well-known wide class can be represented by the 
series 


oe oy G- Ro) y Ge ees ——; n=0,1,2,.... (26) 
k=l 


[or (0) do 
0 


nk 


Thus, the equations (24.1)—(24.3) lead to 


h, GY) 
i. oe hy Gy yi a LS Oe Siena, 
Dee ENCES EAC i ee ONO te a (2734) 
ee Lar Srl 1 Se teed 2... 
wherein 
[a Ryx(o) do g(o, «) de 
1 0 0 
Gh 2m r 
Jo Ri ,(¢) do 
0 
: a (27-2) 
r COSnH & 
(k) fe R,x(9) do | g(a, a) sinn a pe 
Gai ee ee ime ; 
Geiss x a 
“<, / o R2 (0) do 
0 
it Ne ene ts eee | 


Using (21.2), (22), (23), (27.1) and (27.2) gives the wanted solution 
Uy, = U,(0, p, z; t) of our second partial problem 


r 27 


om | Je Ry,(2) do | g(a, a) da 
hy = eG gms ey. : : 3 Ro x(@) Zor (2) 
TT Zoals) + hs Zox(s)[o B3 ,(2) 40 
: + (28) 
; Qn 

= [oe R,, (0) aoe. a) cos” (py — a) da | 
a) a Rx (0) Znx(2) 
py Pie, inf Rial | 


The expressions R,,,(0) and Z,,(2) are defined in VASSAR 


60 VAcLAV VopiéKA ZAMP 


We now proceed to accomplish the operations indicated in (28). Before all, 
we have by (21.2) 


(hy + he) s 6, Cos2 6, , 
s (6,, COS6,;, + hy s sin d, x) 


+ (hy hy s?— 62,) sin2 6,, 


, 


Zin(8) + hy Zy4(S) = | (28.1) 


Bese LOT Sta el nee a 
Now, many technicians will perhaps find useful the following consideration 


making possible to calculate the other integrals entering in (28). If y = y(%) 
is a solution of the equation (x + 0, p > —1) 


KE ap RY. Oe AE ey ee 


then we have to determine i! BW Ags 
Multiplying the above equation by 2 y’ gives 


ae 


pas Pins! Sl a 2 v2 2 
(xy)? = (2 a2 + pa) SY 


dx 


and from this we obtain (K is a constant of integration) 


1 Q ayl2 2 ay2 2 2 ay* Pte 2 
Sey =p y*— 227? Kl =p % oF dx = wtye2fxy um 


Thus, the wanted integral is 


1 
[ey ax = oy [eet x8 + £9) 99 — Py S K, 
Applying these results and taking into account the facts 
R’,() thy R, (7) =9, R,,(0) = (=n 0); bre 2e - n 


we find successively 


Y 


Jo R? (0) do = ee ( ots nt) R? ,(c) — o? R,3,0)| ; 


0 


Me : ‘ 
= Ze, Gat m— her?) Real) 


y2 
ee Y & (En Bee hy 3) He (E,¢) , 
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or finally 


2 
|e R,(0) do =" — (2, + m2 Sh? 1?) P(e, ,); | 
; (28. 2) 


= 1, Zoe = 0, Ta... 


Recapitulating our results the solution ~, of the second partial problem (1), 
(6) is determined by (28)—(28. 2). Herein R,,,(@) and Z,,,(z) are defined in (21.2), 
the values ¢,, denote the roots of the first equation (20) and 6,, are given by 
the formulas 


9 ) s2 
Ean =~ efonte ais eh tee ee ee ee (28. 3) 


§ 4. Determination of u, 


Here we can proceed in two different ways. Firstly, it is obviously possible 
to repeat almost without change the procedure applied at determining the 
partial solution w, and secondly we can obtain the wanted function u,(0, @, z; 2) 
without any calculations directly from (0, @, z;#), ie. from the expression (28), 
writing h(0, y), @2, —h,, —Ah, and —s in place of g(o,@), @,, h,, A, and s. 
This is evident by comparison of the relations (6) and (7). 

Applying the second way gives at once 


8 
‘acces. 
eo. 
Q 
Ee 
So 
a 
aay 
& 
Q 
Peg 
> A 
a 
BS 
X 
R 
— 


"| The Zyx(—s) — Zya(—s)] [0 RB, (0) do 
0 


we Ryrlo) Zual2) / 


Y 2% 5 

’ [eR glo) do [ h(o, a) cosn(p — a) da | | 

"Eh Zyx(—8) — Zin(—s)l [fo R2,(a) do | 
0 


Ss 


En : 
R,,,(@) =e (£28 0), Zn (2) = Sn ~ COS On, + hy S Sindy, : (29.1) 


Re 2 Py on 2 0, 1) 2, o 


E : é 
On p COS On x (1 — =) + hy 5 Sin by x (1 -=) 


and from this 
; (hy + hg) 5 O,~ COS2 O,— + (Ay he 8? — Of,) Sin2 6,4 
hy Znn(—S) Zak s) = e (0,6 Cos ae = hy s sin Ong) , (29.2) 


pales oe gee OOL 20s. > 
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The integrals entering in the denominators of (29) are again given by (28.2), 
e,, denote the roots of the first equation (20) and the values 6,, are defined 
by the formulas 

= — -* V0, + i 95, pala 9 LP Me 0 1,2 


a2 


§ 5. An Important Special Case 


Our solution implies most problems of this kind treated occasionally in 
literature. We expressely notice the case where the surrounding temperature 
fluctuations have constant amplitudes, say Ay, Ay, Ao. 

Here the resulting heat waves in the solid are given by u = Uy + 4, + Ug, 


wherein 
ip IEAGK Is 
ty = hy? Ay eet 5 = : co eta (30) 
[Ro ¥ Lo(&m) + &p Ly (en)] [ ZR(2) at 


= 


hy = ty Ay? DO 1,“ 0) Zyx(2), 31) 
k=1 i 
[Zga(s) + hy Zoq(s)]/ o R3,(0) do | 
0 
A o Ro,(o) do 
thy = hy Ay ei! y ™ 1,(2* 0) Zox(2)- (32) 


iro eeZe. (sie Z,(—s)] [0 R3,(0) do 


The expressions Z,(z) from (30) are defined in (12), if Z3(C) df is given by 


(19.2) and for {[Z4(0) dé one obtains readily 


2 (hy On + he bp COS2 6 hy h 1 

7 dt — 5) Mt oe 2% bt Ay hg s sin2 dx) 

fi w(6) a OF, (6, CoS 6, + hy s sin d,) ; (33) 

The values 6, are positive roots of the equation (10), ¢, are defined in (19.1). 
As for the expressions entering in ot) LZorlZ ¥ are determined in (21.2), 


Zox(S) + Ay Zox(S) in (28.1), the integrals J o Rj,(c) do are given by (28.2) and 
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[ & Ro,(0) do is found, by means of the relations 


" , eR , 
& Ry,(9) + Ro, (0) — a, ea Rox(o) = 0, Rog(r) + ho Rot) a0 


to be equal to 
[o Rj, (0) do = — “a Talfor)i -#=1,2,.... (34) 


d 


The roots ¢), of the equation 
é 1(e}+ hor (ce) = 0 


and the values 69; are connected together by (28.3). 

And finally, Zo,(z) in (32) are defined by (29.1), hy Zyx(—s) — Zy,(—3) by 
(29.2) and 69, can be found from (29.3). All the other expressions and quantities 
entering in (32) are the same as in (31). 

Although no substantial difficulties arise in the effective evaluation of each 
concrete case, belonging calculations are nevertheless very toilsome, on account 
of a comparatively general conception of our beginning problem (1)-(3). 
Fortunately, it often can be assumed in technical practice that the surrounding 
temperature fluctuations have the same constant amplitudes and frequencies 
and that also the heat transfer is the same at all the boundaries of the body in 
question. Denoting those quantities by A, w and h respectively, the above 
general theory leads easily to the following expressions for the resulting 
temperature distribution « = u, + u, within the solid: 


wean iwt Py 5 On 2 
U,=%,(0,2;t) =2Ahre SD TTT is T= ss = 
; ° sin 6 ge o>) 
PD, = D,,(h s) = ine weeds ae 
5, ( T jaa gf sx) 
Z wok Pox ex 
ty = the (0,231) = 2A he ae Rei eas aal. I= 0) cosy 2, 
h 
®,, = ®,, (hr) = Ji(€x) tes a F (36) 


ey JF (ex) + FPen)] (eB + 4? 7?) Soler) 


2 ma) 
B= |—GE + 2) 
the values 6, in (35) and é, in (36) denote positive roots of the equations 


dtgd=hs, efile) —Ar fle) =0. (37) 
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The other expressions in (35) and (36) are to be computed, to the best of 
the author’s knowledge, separately for each concrete technical case. Calculations 
are made slightly easier by the fact that #, contains only even powers of & 
and similarly u, even powers of A,. 


§ 6. Final Remarks 


1. The solution wu of our problem (1)-(3) is given by (4), (19), (28) and (29). 
It is not difficult to recognize the wave character of each component wy. 

2. Naturally, our conclusions hold only for 4) > 0, 4; > 0 and hy > 0. For 
all coefficients # of heat transfer may always be supposed to be positive and - 
the limiting process 4 > 0 generally does not lead to correct results for the 
case, where the respective boundary is impervious to heat. 

3. Our deductions can evidently be applied to several important industrial 
questions. 


Zusammenfassung 


Vorliegende Arbeit bringt eine vollstandige theoretische Berechnung des 
Warmeregimes im endlichen Kreiszylinder, dessen Oberflache periodischen 
Temperaturschwankungen mit lokal veranderlicher Amplitude ausgesetzt ist. Die 
Ergebnisse sind auf eine Reihe technisch wichtiger Falle anwendbar. 


(Received: May 8, 1956.) 
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Knicken eines gestreckten Gelenkstabzuges 


Von Hans-H. Gasser, Lungern (Schweiz) 


Wird ein Stab zentrisch gedriickt, so knickt er unter einer bestimmten kriti- 
schen Last aus. Dabei kann die Knicklinie einen Wendepunkt haben, was bedeu- 
tet, dass an dieser Stelle wahrend des Knickens kein Biegemoment entsteht. Man 
kann sich daher fragen, ob dort ein Gelenk angeordnet werden darf, ohne dass 
dadurch die Stabilitat beeintrachtigt wird. / 

Die folgende Untersuchung zeigt, dass das in einem Fall tatsaéchlich méglich 
ist, in andern Fallen aber nicht. 

Zu diesem Zweck soll zunachst das Knicken eines Systems nach Figur 1 
untersucht werden. Hier sind zwei Falle zu unterscheiden: a) Der Stab 1 knicke 
vor dem Stab 2 aus. Fiir 1 liegt somit der einfachste Euler-Fall vor, und da die 
Belastung von 2 auch nach dem Knicken zentrisch bleibt, wird das Gelenk 
zwischen 1 und 2 nicht verschoben. b) Der Stab 1'sei so kraftig, dass der Stab 2 


unter geniigender Belastung ausknickt. Fiir diesen Fall soll nun die kritische 
Last berechnet werden. 
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Der Stab 1 habe die Lange a, der Stab 2 sei homogen und prismatisch, unten 
total eingespannt und habe die Lange /. Die Differentialgleichung der elastischen 
Linie kann nach der Energiemethode ermittelt werden. Diese beniitzt die Tat- 
sache, dass fiir das Gleichgewicht des Systems die Summe der potentiellen 
Energie der Belastung und der Deformationsenergie stationar ist. 

Die Verkiirzung der Stabe unter der Belastung beeinflusst das Knicken in 


der zu beobachtenden Naherung nicht. Demnach ist die potentielle Energie der 
Belastung 


1 
V =—Pla— a? — y2(l)]—P , : VL y’2(x) | 
iar 


oder bei Beschrankung auf kleine Ausbiegung 


1 
phot fe, | 
vV=—plee so eae aa 
0 


Die Deformationsenergie ist 


Die Variation 6{[V + U]= 0 ergibt die Differentialgleichung 


Py" + EJ yl 0 Figur 1 


2 A Mechanisches System, 
sowie die dynamischen Randbedingungen 7 


7 (=O (Oram Gelenk) 


4 [ yl) =p v'(0)| + EJy"(l)=0 (Gleichgewicht der Querkraft im Gelenk). 
a 


Nimmt man noch die zwei kinematischen Randbedingungen 
y(0)= 0 und y’(0)=0 (totale Einspannung unten) 


dazu, so kann die Differentialgleichung 4. Ordnung gelést werden. Ihr allge- 
meines Integral ist 
y(*) = Acosxx+ Bsinu. 7+ Cx+D, 


wobei x? = P/E J ist. Setzt man y(v) in die 4 Randbedingungen ein, so erhalt 
man ein 4gliedriges, homogenes Gleichungssystem fiir die Unbekannten A, B, 
C und D. ' ; 

Soll eine nichttriviale Knicklinie méglich sein, so muss die Determinante des 
Gleichungssystems 


1 0 0 ‘st 

| % 1 0}. 

Cd aa | cosxJ sinxl 0 0 | 
| 0 0 eee ah 


verschwinden. A(x) = 0 fiihrt auf die transzendente Gleichung 


“(a+ /l)=tgxl. 


ZAMP VIII/5 
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Mit dem Verhaltnis (a + /)/7 = mw folgt 


wul=texl.| (I) 


Aus Gleichung (I) kann x bestimmt werden, womit die kritische Last 


| Pari, = ES | (11) 


berechnet werden kann. 
Die Gleichung (I) kann nach Figur 2 graphisch gelost werden. 


Figur 2 


Graphische Lésung der Gleichung (1). 


Ballin co a) ] == (055 (a 0 oe) sO aen (ge Oe P beliebig ; 

b) a>o; mle 9t) 2: Euler-Fall; 
HalliZc = cols | waaid 10); i) 2 oO; 
alles sei —eles a 0r x= 4,49 (P, = 20,16 Bf /l?): Euler-Fall; 
Pallas LS 7 0; Sa — 7; 4,49>xl>a7: Stab 1 ist Zugstab; 
all 3 pies Os Ci. son Euler-Fall; 
Ballon <= 0: (ee ahh GS “lS 7/2. 


Vergleiche die entsprechenden Knicklinien in Figur 3. 

In einem zweiten Schritt soll das erweiterte System von Figur 4 untersucht 
werden. Weitere Gelenke kann man nicht mehr einfiihren, da sonst das System 
fiir beliebige Werte von P unstabil wird. Auch hier bleiben die Gelenke in Ruhe, 
wenn man annimmt, der Mittelstab knicke zuerst aus, da fiir diesen der einfachste 
Euler-Fall vorliegt. Interessant ist wieder der zweite Fall, wo die Stabe 1 und 2 
ausknicken. Beide seien homogen und prismatisch. Sie haben die Langen /, bzw. ~ 
/, und die Biegesteifigkeiten (EJ), bzw. (EJ),. Es ist noch zu bemerken, dass das 
Ausknicken der Stabe nur in einer gemeinsamen Ebene médglich sein soll, da 
sonst Torsionsmomente auftreten wiirden. 

Die Differentialgleichung kann analog gefunden werden. Das Variations- 
problem liefert hier die beiden Differentialgleichungen 


P yj (4%) + (EJ yi (4%) = 0, 
P yy (4%) + (EJ)o 2" (%2) = 0 
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Fall: 1b 


Figur 3 Figur 4 
Zur Diskussion der verschiedenen Knickfalle. Das erweiterte System. 
und die dynamischen Randbedingungen 
Villy) = 9; yo(lz) = 0; 


(s gee 4 yi(h)| + (EJ), ¥'() = 9, 


P [2a — lh oF va(l)| + (EJ)2y2'(22) = 0. 


Dazu kommen noch die kinematischen Randbedingungen 
y,(0) = y1(0) = y,(0) = y3(0) = 0 (totale Einspannung). 


Mit den Abkiirzungen P/(E J), = x? und P/(E J), = § sind die allgemeinen 
Integrale 


Wy, = A, cosx, ¥; +B, sin, 4, + Cy*#,+ D,, 
Yo = A, cOSH#,%, + By sinx,%,+ Cy%, + Dz. 
Die 3. und 4. Gleichung der dynamischen Randbedingungen lassen sich nach 
Einsetzen der Ausdriicke fiir y, und y, auf 
C, = C, (= C) 


und 
C (ly + lp + a) + Di + D, = 0 


umformen. Die zweite Gleichung bildet mit den 6 noch nicht verwendeten Rand- 
bedingungen ein 7gliedriges, homogenes Gleichungssystem fiir die 7 Unbekannten 
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A,, Ay, B,, By, C, D, und D,. Die Determinante 


1 0 0 0 0 iby 
0 1 0) 0 0 0 tf 
0 0 Hy 0 1 0 0 
A= 0 0 0 Ho 1 O04 
cos x, 1, 0 sin x, J, 0 0 Oy | 
0 COS Hy 1, 0 Sin x» 1, 0 O 0.) 
0 0 0 0 (i Ly. a ee 1 
gleich Null gesetzt, ergibt die Gleichung 
1 1 
t, + ily + a= ape ke = yrs Sor pee 
Mit J, + J, + a = L (Gesamtlange), L/J, = w, und L/l, = pe folgt 
1 tgmyh4 , 1 tg x5 1, ie | (IIT) 


es xy dy Ie Hy ly 


Das Verhaltnis der Winkel x, /,/x. /, ist durch das System gegeben, da 


WPIEDs. 
2 VPKME De 


ay 1, ry l, eee 
: He ly ly (EJ\i 
ist. 


Fiir die graphische Lésung der Gleichung (III) zeichnet man am besten die 
beiden durch mw, bzw. my dividierten Funktionen tgx, J,/x,1, und tg x, 1,/x% 1, 
libereinander und sucht die Punkte, an denen die Ordinatensumme gleich 1 ist. 

Die Funktion tgx//x1 ist zum besseren Versténdnis der folgenden Diskus- 
sion in Figur 5 dargestellt. 


und somit 


Figur 5 
Zur graphischen Lésung der Gleichung (III). 
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Ist a positiv, das heisst, ist der Mittelstab auf Druck beansprucht, dani ist 


a ite 


und die Lésungen x, /, und x, /, liegen beide im Intervall 0< xl< a/2 (Figur 6, 


|e 


Fall 1). 


Figur 6 


Die verschiedenen Knickfalle. 


Ist “, = “, und x, = x, (symmetrisches System), so liegt die kleinste Lésung 
xl ebenfalls im Intervall 0< x/< a/2. Die zu ihr gehdrende Knicklinie ist 
antimetrisch (Fall 2). Die nachsthGhere Lésung ist x/ = 2/2. Der Mittelstab 
bleibt senkrecht, und fiir die Stabe 1 und 2 liegt ein Euler-Fall vor (Fall 3). 

Erst fiir a = 0, wo 


ist, ist die kleinste nichttriviale Lésung 7 = a/2 (Fall 4). 

Liegt das nunmehr einzige Gelenk nicht in der Mitte, dann ist x/ des einen 
Stabes knapp iiber, jenes des andern Stabes knapp unter 2/2. Die genaue Form, 
insbesondere das Auftreten eines Wendepunktes kann erst anhand der ausge- 
werteten Funktionen y, und y, ermittelt werden. 

Ebenso ist fiir a< 0x,/, und x, /, groésser bzw. kleiner als 2/2 und bei Sym- 
metrie gleich 2/2 (Fall 5). ' 

Mit den Gleichungen (I) und (III) kann nun die am Anfang gestellte Frage 
beantwortet werden. 

Fiir den Fall, wo das eine Ende der Stiitze total eingespannt und das andere 
Ende gelenkig gelagert ist, betragt bekanntlich 


EJ 
ay 


Prarit. = 20,16 


und der Momentennullpunkt ist 0,30 L vom eingespannten Ende entfernt. Ordnet 
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man dort ein Gelenk an, dann ist in Gleichung (1) 


a+l Je 


a = eee 
be j 0,30 L : 


einzusetzen. Sie ergibt als kleinste Losung x / = 1,35. Hieraus nach (II) 


352 Sr 


f oe (0,30 L)2 


krit. 


EJ 
=x? EJ =- E Jie 20,16 = 
Der Knickvorgang wird hier tatsachlich durch Einfiihren eines Gelenkes im 
Momentennullpunkt nicht veradndert. 
Anders liegen die Verhaltnisse bei der beidseitig eingespannten Stiitze. Ihre 
kritische Belastung ist nach EULER 


4n° EJ 


Pi L2 


Fiihrt man a) nur in einem Wendepunkt (hier im Viertelspunkt) der Knicklinie 
ein Gelenk ein, dann ist /, = L/4, 1, = 3 L/4 und a = 0. Somit pw, = 4, My, = 4/3 
und x, J,/x5 1, = 1/3. Gleichung (IIT) lautet 


1 tgm 3 tg3%,1, 


. = i 
4 Ze thy 4 3 x, 1, 
oder 
tgx,1,+ tg3x,14,=4%,1,. 
Eine Lésung ist x, /, = 2/2. Daraus folgt 
Pee 2G mu ET 4 3? 
eho) we dE z: ese bls pat 
Pai IVES |s urge epee Tapes: 


in Ubereinstimmung mit obigem Wert. Es gibt aber noch eine kleinere Lésung 
im Intervall 2/3 > %/> a/6. Der Knickvorgang wird in diesem Fall verandert, 
da das System schon unter der kleineren kritischen Last ausknicken wird. Bei 
Variation der Gelenklage kann die Knicklast erhdht werden und erreicht den 
Héchstwert, wenn das Gelenk an einem Ende liegt, was dem vorgangig behan- 
delten Fall entspricht und daher hier sinnlos wird. 

Wird b) in beiden Wendepunkten des beidseitig eingespannten Knickstabes 
ein Gelenk angeordnet, so ist 1, = 1, = L/4, a= L/2 und somit w, = uw, = 4 und 
#, 1,/%, 1, = 1. Gleichung (III) lautet 

Lite ccna Deatearn Us 


at nity, eres 


=1. 


_ Die Losung x, 1, = x, /, = 2/2 stellt wiederum den Euler-Fall dar. Aber auch 
sie ist nicht die kleinste Lésung. Das Verhaltnis « = 4 braucht demnach nicht 
das giinstigste zu sein. 

; Zum Schluss soll fiir den zuletzt beschriebenen Fall jenes ~ gefunden werden, 
bei dem P,,,, maximal wird. In diesem Idealfall wird die kritische Last fiir den 
Mittelstab und das ganze System gleichzeitig erreicht. Der Mittelstab, der die 
Lange 

a=L—h~h=L—-2==2(1- “ 
Mu i 
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hat, knickt nach EuLER unter 


Pp dak ot Rall ney 


krit. = fies 


~ LYS of yy? 
aus, das ganze System aber nach “Gleichung (II) unter Pi = «2 E J. Durch 
Gleichsetzung erhalt man x = a/L (1 —2/u), woraus x1 = “ L/M = n/(u — 2) 
oder uw = (a + 2x 1)/x1 folgt. Gleichung (III) heisst etwas umgeformt 


mu eingesetzt ergibt 
Z 
mt Sule ee 
woraus x / = 1,23 als kleinste Lésung berechnet wird. Daraus folgt u = 4,55, 
und die Lange des Mittelstabes wird 0,561 L. 
Die kritische Belastung ist 
1,23? 31,3EJ] 
Or ee 


Prrit. =x FE] 


somit im giinstigsten Fall noch 79% von derjenigen der durchlaufenden Stiitze. 


Summary 


The buckling of a multi-pinned bar (see Figures 1 and 4) has been investigated. 
In particular, the influence of pins along the bar at the zero-moment points is 
considered. 


(Eingegangen: 30. Juli 1956.) 


The Second Fundamental Problem of Elasticity Applied 
to a Plane Circular Ring 


By Wittram A. Gross, Murray Hill, New Jersey, USA?) 


Introduction 


H. REISssNER?) has discussed the problem of a plane elastic ring under dif- 
ferent boundary conditions. One of the terms in his general solution appears 
to be in error. 

The general solution of the second fundamental problem, that in which 
boundary displacements are specified, for a plane circular ring is given here. 
The complex variable method of N. I. MUsKHELISHVILI') is used. The solution 
of a particular problem, that of a ring with fixed outer boundaries and a rigidly 
translated inner boundary, is compared with the solution of H. REISSNER. 


1) Murray Hill Laboratory, Bell Telephone Laboratories, Inc. 
2) H. ReIssNnErR, Jb. wiss. Gesellsch. Luftfahrt 1928, 126 fies 
3) N. I. MuskHexisuvit1, Some Basic Problems of the Mathematical Theory of Electricity, 


Radok translation (P. Noordhoff Ltd., 1953), p. 218 ff. 
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Analysis 


It may be demonstrated [MUSKHELISHVILI%), p. 138], that, if g(z) and yw(z) 
are holomorphic in a plane region G, admissible stresses and displacements in 


polar coordinates are given by 


0,, + O99 = 4 Rey’, (1) 

o,,—i6,,=2Req— 2 Fo" + ¥), mee 

2 (u,+ iu) =e [xp— zy — ¥), (3) 

in which z= ve'", wis the modulus of rigidity, and x = (3 — v)/(1 + v) for plane 


stress or 3— 4v for plane strain, v being Porsson’s ratio. 
Let the inner radius of the ring be « and the outer radius £. Represent .the 
known boundary displacements by the Fourier series, 


21 (uy im) 3) Ane, (4) 
onv =a, and oe 
2 uw (u, + iu) = S' Be”, (5) 
ony = 6. Set ~ 
g =A logz+ )" a, 2”, (6) 
—00 
and Me 
y= —x A logz+ 3"d, 2". (7) 
co 


The first terms in equations (6) and (7) imply single valued displacements 
everywhere in the ring. The coefficients are determined from equations (3), (4), 
and (5). j 
A Bie Be Ae (Bpk Ae (8) 
2x8 (x2 +f) log® — 2 (8% — a 


2% A log ® + oe ee 


ag = 2 (p2 — 0.2) ) (9) 
b_» = a2 [xa,a2— A— A], (10) . 
% diy —' by = A_,— 2% A loge + 20°4,,' (11) 


and 


(a= 2) (B*— 24) (a-" Ay y— BP" By») + (a Ag _g— BM By.) x (a-™— ft-am) 


22 (q4-2n B4-2n) (2% — £2”) +n (n— 2) (B2— a?) 2 ? 


b= CAL Ae xy Cie (2> n) ay» ae ile (13) 


for all different from 0 and 2. a, and b, do not affect stresses and enter into 
displacements only in a combined relation like that of equation (11). 
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Example 


Let the inner boundary of a plane circular ring be rigidly displaced 6 in’ the 
direction 9 = 0, Fix the outer boundary at the unloaded position. Hence on v = B 


2 Ul (UperbtUe) = 0, 
and, for 6 small compared to «, 
2p (u, + iu) =2 nde. 


It follows that 4_,= 2 w 6 and all other 4, and B,, are zero. Only the following 
real coefficients are different from zero: 


ct Ballitin. Ox mp (a? + B?) 


~ arate onde: (14) 

ma (one te, RE log—— s((6" — at) 
A : 
72 (a? + BB) C3 

—a2 B24 : 
eS eee 16 
xa — b= 2.4 |x log p : & (17) 
pele x (a? + B) |. 
It follows from equations (2), (3), and (4) that 

6,, = [(3 + *) A va + 2 ag ¥ + ea le cos, (18) 
Og = [(1— *) Ar 1+ 6a,7— 26,7 *) cos6, (19) 
65=[(1— x) Ar 1+ 2a,r+ 2b_or *Jsin@, (20) 


u, = (2 w)*[2%A logy + (x — 2) ag 7° + xa) — by - A—bd_gy*] 0086, (21) 
Uy = (2 w) 1 [—-2% A logr + (x + 2) a, 7 — xa) + by — A —b_g7"] sind. (22) 


For the same problem (REISSNER®), p. 130), gives, in the present notation, the 
complex displacement 3 
log — 

Y 
B ’ 


o4 


i0 | 


U, + tg = 0e- (23) 


log 


which obviously satisfies the boundary conditions, but not the equilibrium 
equations. 


_ Zusammenfassung 


'Ein ebener Kreisring mit vorgegebenen Randverschiebungen wird mittels 
der komplexen Methode nach N. I. MUSKHELISHVILI behandelt. Es wird ein ein- 
faches Beispiel gelést und die Lésung mit der fehlerhaften von H. REISSNER 


verglichen. 


(Received: July 31, 1956.) 
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A Graphical Method for Determining the Coefficient of Viscosity 
of Newtonian Liquids Using an Oscillating Cylinder Viscometer 
II. 

By Att ABDEL Keri [prAuiIM, Alexandria, Egypt’) 


Consider a Newtonian liquid of coefficient of viscosity 7 in g cm~*s~* and 
density @ in g cm~* is placed between the two cylinders. By varying the frequency 
n in c/s of the outer cylinder, the angle g, i. e. the phase angle by which the inner 
cylinder lags the outer cylinder is changed and according to my equation [1]?). 


n= 2ude 24 (4 dg) cot, (1) 
where wen . 
b—a wo2—T +a 
= Se Th . 
3 QO EN Bee ae, at aia 2 


b is the radius in cm of the outer cylinder, a the radius in cm of the inner cylinder, 
I the moment of inertia in g cm? of the inner cylinder, zt the torsion constant of 
the wire in gcm? s-2 and L in cm the depth of immersion of the suspended 


cylinder. 
If tis very small compared with I w?, we have, 
iL 
—3d/R 
n= 2ode : lasrer + 4 e] cote. (2) 


If the end effect ¢ is considered, equation (2) becomes 


Lo pate 
22a*(L + é) 


n= 2ade 24h | +49] cot. (3) 


This equation represents the coefficients of viscosity 7 of non-elastic liquids. 
It therefore follows that if tan m be represented by Y and w be represented by 
X, then the relation between Y and X is linear and may then be written. 


Y= Xx, (4) 


where 


I 
= -34/R 
Reraae Bae 


& depends on the dimensions of the apparatus used and the density of the - 
liquid under test. 

It therefore follows that the slope of this line is K/n. The figure shows Mar- 
KOVITz et al, [2] results taken with standard viscosity oil obtained from the Na- 
tional Bureau of Standards (OB-5 oil). 

The obtained graph is considered to be satisfactory. The agreement indicates 


that there is no end effect; e = 0, and therefore the effective depth of immersion 
is the actual depth. 


1) Physics Department, Faculty of Science, University of Alexandria. 
2) Numbers in brackets refer to References, page 75. 
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The following table shows the results for the liquid examined: 


7 in Poises 
Slope K/y~—/ |-—-—_-——_-—___ 
| from slope | from a paper 


0-470 x 10-2 85-1 | 91-6 


It is seen from the figure and table that equation (3) gives a representation 
of the experimental facts. 


Also the figure shows that » approaches 2/2 as the frequency is increased. 


tan Y 


2 4 6 


The variation of tangy w for National Bureau of Standards oil OB-5. 


REFERENCES 


[1] Arr A. K. Iprauim, ZAMP (in preparation). 
[2] MarKow17z et al., Rev. Sci. Inst. 23, No. 8 (August 1952). 


Zusammenfassung 


Es wird ein zweites graphisches Verfahren dargestellt, welches gestattet, die 
Viskositat ziher Fliissigkeiten zu ermitteln. 


(Received: June 4, 1956.) 
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On the Motion of a Projectile in the Atmosphere 


By V1-CuHenG Liv, Ann Arbor, Michigan, USA?) 


Abstract 


The equation of rectilinear motion of a projectile which moves in an atmos- 
phere, of which the density decreases exponentially with the altitude, is solved. 
It is found that. the velocity of the projectile can be expressed explicitly in terms 
of confluent hypergeometric functions. This theory is applied to treat two 
specific problems: (1) the flight analysis of a sounding rocket during the free- 
flight period and (2) the calculation of ambient temperature from the trajectory 
of a spherical projectile. 


1. Introduction 


The theory of flight of a ballistic projectile in a resisting medium dates back 
to Newton. In the meantime, various approximations for the flight parameters 
which are involved in the equation of motion have been used [1]?). The most 
important ones are the aerodynamic-drag coefficient of the projectile and the 
ambient-air density. The former is usually considered as functions of Mach 
number, the latter as functions of the altitude. 

The theory of flight of a sounding rocket during the free-flight period deals 
essentially with the same problem as that of a ballistic projectile. The only 
difference is that the air-density factor is more critical because of the extreme 
altitude attained by the former; hence, a more realistic approximation to the 
ambient-air density is necessary. In view of the multiplicity of parameters in- 
volved, the equation of motion of the projectile (or rocket) is generally treated 
either by a laborious method of step-wise integration or by the analog computer 
technique. 

A recent experiment [2], in which a falling sphere is used as an ambient- 
temperature probe in the upper-atmosphere measurements, stimulates new 
interest to the trajectory problem. Measured trajectory data of a spherical pro- 
jectile are used to determine the ambient-air density, from which the Lespugiash: -air 
temperature can be derived. 

In this note, new approximations to the aerodynamic-drag coefficient and 
ambient-air density, which are believed to be more accurate than previous repre- 
sentations [1], are introduced. Particle drag, which is caused by the presence of 
aerosols in the atmosphere, is allowed in the present analysis. This, of course, 
warrants consideration only in some special cases. 

It is found that the projectile velocity is expressible explicitly in terms of 
confluent hypergeometric functions. The significance of this result, which is 
believed to be new, is twofold: (1) In treating ‘direct’ problems, for instance, 
like the trajectory analysis of a sounding rocket, the present theory can be used 
to save the labor of step-wise integration if tabulation of the particular confluent 
hypergeometric function is available. (2) In the case of an ‘indirect’ problem 
such as the computation of the ambient temperature from a falling-sphere 
trajectory, it can be applied with advantage because of the elimination of the 
numerical differentiation process which is involved in the original theory [2]. 

1) High Altitude Engineering Laboratory, University of Michigan. 

2) Numbers in brackets refer to References, page 81. 
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2. Significant Forces 


The word ‘projectile’ is used here referring to any axially symmetric body 
in free flight. We restrict our discussion to the case with negligible spin and yaw. 

The force system acting on a projectile includes: (1) the gravity mg; (2) the 
buoyancy @g V, where V represents the volume displaced by the projectile, and 
g the ambient-air density; (3) the aerodynamic drag D; (4) particle drag d, which 
is caused by the presence of aerosols in the atmosphere; and (5) the inertial force 
(m+ Am) dv/dt, where Am, Munk’s apparent additional mass, is equal to one 
half the mass of air displaced by a spherical projectile; and du/dt, the projectile 
acceleration. 

The buoyancy and the inertial force of the apparent additional mass are 
significant only for projectiles of the balloon type. Though, under sea-level 
conditions, the particle drag due to water droplets, for instance, may be insigni- 
ficant when compared with the aerodynamic drag, it is conceivable, however, 
that in the cumulus clouds at high altitudes, the drag due to water droplets may 
become significant. While the water-droplet content depends primarily on the 
ambient temperature, the ambient-air density decreases with the altitude expo- 
nentially. 

It is not intended here for us to become involved in detailed computation of 
the particle drag. Rather, we shall investigate its effect on the projectile motion 
since the particle-drag term in the equation of motion follows a functional 
dependence upon the altitude, which is, in general, different from that of the 
aerodynamic drag. As a first approximation, the particle drag is prescribed as 
d= C,0A v?/2, according to the Newtonian concept of drag. The particles are 
assumed to be mass points at rest; hence, Cz should depend only on the position 
and shape of the projectile and be independent of its size A und velocity v. The 
particle concentration is represented by o. 


3. Approximate Representation of the Aerodynamic-Drag Coefficient 


It can be demonstrated with a dimensional analysis that the aerodynamic- 
drag coefficient, defined as Cp = D/(g V?-A/2), for geometrically similar projec- 
tiles depends primarily upon the Mach number and Reynolds number. For 
projectiles moving at high speeds, the dependence of Cp on the Reynolds number 
is much less significant than that on the Mach number except in cases for which 
the Reynolds numbers are low, such as those at the high altitudes. In those 
cases, however, the aerodynamic-drag force is usually negligible as compared to 
the gravity. 

For a specific projectile, the variation of Cp with Mach number M, in general, 
shows distinctly different characteristics in the three Mach number ranges: 
(1) subsonic, (2) transonic, and (3) supersonic. Qualitatively speaking, Cp(M) 
starts with roughly constant values (for different projectiles) then increases 
gradually as transonic range is approached. During the transonic range it in- 
creases until the supersonic range begins. Cp(M) in the supersonic range decreases 
_ with a gradient which is very large at first and gradually diminishes to a small 
value. i 
Reliable values of Cp still rely mostly on measurements from model tests. 
_ Theoretical consideration may, however, serve as a guide to the choice of empiri- 


cal functions for prescribing Cp. 
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Consider separately the compressibility effect in the subsonic flow for slender 
and round bodies. In the former case, to account for the compressibility effect on 
the pressure coefficient, either the Prandtl-Glauert factor (1— M?)—1? or a 
modified version of it (for instance, the Karman-Tsien formula or GOETHERT’S 
generalized Prandtl-Glauert formula) can be used. These factors, when expanded 
in a power series of M, can be expressed in the form 


Cy 4 Cpl tape (3.1) 


after neglecting higher-order terms. In the latter case, as an index of the compressi- 
bility on the drag coefficient of a round body, the impact-pressure coefficient 
at the stagnation point can be shown to be 


M2 
Sty pee (3.2) 


It has also been shown by measurements that the compressibility effect on the 
skin-friction coefficient and the coefficient of drag due to interference between 
body and fins can be considered as proportional to the powers of the Prandtl- 
Glauert factor. 

A large fraction of the drag force on a projectile in the supersonic flow are 
attributed to the frontal shock wave and the base pressure. For a typical body 
like a cylinder with a nose cone, the drag coefficient can be shown to be approxi- 


mately of the form 
C,+ C, M-?. (333) 


The validity of the above expression is expected to extend to the hypersonic 
range [3]. Relatively less is known concerning the theory of drag in a transonic 
flow, which is a mixture of subsonic and supersonic flows. It is therefore sug- 


gested to use 
Cp = C,+ C, M2+.C, M-* (3. 4) 


as a general expression for approximating the drag coefficient of a projectile in 
the following analysis. Cy, C,, and C, in (3.4) are empirical constants to be 
prescribed for each Mach-number range. 


4. Equation of Rectilinear Motion of a Projectile 


Consider a nonspinning projectile which ascends at zero angle of yaw in a 
stationary atmosphere. By collecting terms of forces which have been discussed 
in Section 2, we obtain the equation governing the rectilinear motion of a 
projectile along the vertical axis y. 


1 
am (Co + C, M2+ C, M-?) A ov? 
1 a dv | 
+z Cad ov —egV+tmgst (m+ ils fils gS 
For an isothermal atmosphere with temperature T in altitude interval (y — ‘viele 
the ambient-air density can be expressed as 9 = 0) exp [— (y — y,)/H] where H, 


the scale height, is defined as kT/m, g (k = BoLtTzMann’s constant; m, = mean 
molecular mass of air). 
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To simplify equation (4.1), we introduce 


Z = M*= (yg H)=' 02, * = exp— 270, and 9= 0%, (4. 2) 
where y is the ratio of specific heats of air, and obtain 
dZ B , € 
Ge = 2024+ (Coat £)Z+ (Cra- 54 4), (4.3) 


where 


a= A eH (m+ Am)-', EN RT | As 


60=20,V yy! (m+ Am), e=2my-1 (m+ Am). 


Equation (4.3) is nonlinear and of the Riccati type. It can be transformed 
through use of the substitution 


u = exp [- Je, Zar] (4.5) 
into a linear equation 
d*u B\ du | rene 
age ~ (2 Cot 5) Ge + eC, (a6,- 34+ =) uno. (4. 6) 


Equation (4.6), which has a regular singularity at += 0 and an irregular 
singularity at x = oo, can be converted into the canonical form 


d*p dp = 
eS + b-) SP —ag=0 (4.7) 
with 
RAL OLS ft eS et (6— « C;) = U ex Ome Cyle 
= , Nees oa 2 Seite 2 pina Fo 
Note (4. 8) 
eng a= = (540 p—-aC,e—fa)ar. 


One of the fundamental solutions of (4.7) is the confluent hypergeometric 
function [4] 


(a, b, &) = (a, 0, 8) = Yas, (4.9) 


where 


The other one is 

&-> F(it+ta— b,2— 5; &) (4.10) 
provided } is not an integer. The general solution of (4. 7) with nonintegral b may 
be written as 


® = B, F(a, b; &) +C &-*,A(1+ a—b,2-— 6; €), (4.11) 
where B and C are arbitrary constants. Confluent hypergeometric functions have 


been fairly extensively tabulated ([5j, [6], [7]). 
In view of (4.2), (4.5), (4.8), and (4.11), we obtain 


EE eee 1 =|: (4.12) 
a Gs 


y2 


2 EDswa a “ede 
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Because of the term (d@/d£)/®, one of the arbitrary constants of ® as given in 
equation (4.11) will be eliminated in equation (4.12) for v?. Since equation (4. 3) 
is of the first order, there is only one constant of integration. This is consistent 
with the physical aspect of the problem. The constant of integration is to be 
determined by the use of the boundary condition of the projectile trajectory. 

Consider the restriction on b in (4.11). Since b = — f, where f is defined as 
C,AoH (m+ Am)-1, which, in general, is less than unity, the limitation on the 
values of b is not significant as far as the present analysis is concerned, The 
only notable exception is the case b = 0, which corresponds to the problem when 
the particle drag is null. 

Equation (4.7) with 6 = 0 has been studied by BaTeMAN [8]. The solution of 
(4. 7) for this particular case is known as BaTEeMAN’s k-function. The general 
solution of equation (4.7), when b = 0, is bounded in the neighborhood of § = 0%). 

It is to be noted that when C, = 0, equation (4.3) becomes a linear equation 
of the first order, which can be integrated immediately in terms of elementary 
functions [9]. 

It is also to be noted that the equation of motion of a projectile with the 
assumption Cp ~ M~! (a satisfactory approximation for the wave drag of airfoil 
sections in hypersonic flows) and B= 0 is again a linear equation of the first 
order if v and ¢ are used as a dependent and an independent variable, respec- 
tively [10]. 

It is interesting to note that the function » in equation (4.7) can be expressed 
in terms of Bessel functions provided that the condition 2 C, e = C, # is fulfilled. 
Other special cases of the confluent hypergeometric function will be obtained 
for the solution m when values of the parameters 6 and a in equation (4.7) fulfill 
appropriate restrictions in each individual case. 


5. Engineering Applications 
The result of the present analysis may be applied to the following problems: 


5.1. Flight Analysis of a Sounding Rocket 


A rocket which carries sounding instruments to high altitudes for measuring 
the upper atmosphere is called a sounding rocket. It is usually launched verti- 
cally upward to gain the utmost peak altitude. In view of the simplicity of the 


designed trajectory, it is generally equipped with the least guidance. As a first — 
approximation, it is plausible to assume that a sounding rocket with normal 


performance moves rectilinearly at zero angle of yaw. Hence, the present analy- 
sis can be applied in the free-flight period. 

For a solid-propellant rocket, which usually has large jet-propulsive thrust 
of a short duration, we can estimate the burnout velocity and altitude by neg- 


lecting the aerodynamic drag in comparison to the forward thrust. The classical. 


method is to integrate the equation of motion of a rocket with mass m as a 
function of time ¢: j 


dm dv 
ca +m (ZF +e) == 05 (5.1) 


where C is the constant exhaust velocity. This leads to 


My - 
v= — Cin & (& — ty), (5. 2)% 


0 


3) This is pointed out to the author by Professor Ropert C. F. BARTELS. 
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2 m 1 
Vy mes (in) dt — erty — ty)? (5.3) 


where subscripts 0 and b denote initialand burnout condition, respectively. The 
above calculations can be improved through an iterative process in regard to the 
aerodynamic-drag correction. 


5.2. Ambient Temperature Computation from Trajectory Data of a Falling Sphere 


The idea of using a free-falling sphere, of which the trajectory can be accu- 
rately measured, to probe the upper atmosphere for the purpose of measuring 
the ambient temperature was originally conceived by Jones[2]. This was 
carried out by incorporating a miniature Doppler receiver-transmitter in an 
inflated sphere, which was released from a sounding rocket near its peak altitude. 
The sphere trajectory is determined by a DOVAP system [11], which consists of 
a transmitter and an array of receivers on the ground. 

In an alternative scheme, as proposed by the author, the ground stations of 
transmitters and receivers are replaced by.a mobile station which has a trans- 
mitter and a receiver just like those used in the DOVAP system except in 
miniature sizes. The mobile station is incorporated in a body, with negligible 
drag-to-weight ratio, which is to be released simultaneously with the sphere from 
the rocket. With this scheme, the velocity of the sphere, relative to an apparently 
null-drag trajectory or pseudo-vacuum trajectory, is obtained from Doppler 
cycles which are to be telemetered from the mobile station and recorded along a 
time scale. It is assumed that both the mobile station and the sphere fall verti- 
cally with no relative tumble. 

One serious drawback of the falling sphere method is the need of numerical 
differentiation of the velocity data in the process of computing ambient-air 
density according to the original theory [2]. As it is known, numerical differen- 
tiation is a notoriously inaccurate process because it exaggerates the irregulari- 
ties of the numerical function whose derivatives are to be determined. This 
process can be avoided if the present theory of projectile trajectory is applied in 
step-wise calculation of the scale height H, which is proportional to the ambient- 
air temperature T. 
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Zusammenfassung 


Die Bewegungsgleichung eines vertikal ansteigenden Geschosses wird gelost 
fiir den Fall, dass die Dichte der Atmosphare exponentiell mit der Héhe abnimmt. 
Die Geschwindigkeit des Projektils kann dabei in konfluenten hypergeometri- 
schen Funktionen ausgedriickt werden. Die Theorie wird hierauf auf zwei spezi- 
fische Probleme angewendet: 1. Fluganalyse einer Vertikalrakete wahrend der 
Zeit des schubfreien Fluges und 2. Berechnung der Umgebungstemperatur der 
Bahn eines spharischen Projektils. 
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Numerische Behandlung von Differentialgleichungen. Von L. CoLLaTz. 
Zweite, neubearbeitete Auflage (Springer-Verlag, Berlin 1955 [Grundlehren der 
mathematischen Wissenschaften, Band 60]). 526 S., 118 Abb.; DM 56.—/59.60. 

Die zweite Auflage des Handbuches iiber numerische Lésungsmethoden fiir 
Differentialgleichungen unterscheidet sich von der ersten Auflage vor allem 
durch Aufnahme der in der Zwischenzeit gewonnenen neuen Erkenntnisse, ins- 
besondere der Instabilitatsphanomene bei numerischer Integration. Ausserdem 
wurde die Einteilung des Werkes etwas verandert, indem alle mathematischen 
Grundlagen im 1. Kapitel zusammengefasst sind. Es folgen wie friiher: II. An- 
fangswertprobleme; III. Randwertprobleme bei gewéhnlichen Differentialglei- 
chungen; IV. Anfangs- und Anfangsrandwertaufgaben bei partiellen Differential- 
gleichungen; V. Elliptische partielle Differentialgleichungen; VI. Integral- und 
Funktionalgleichungen. A. Rutishauser 


Mécanique vibratoire. Par Rospert Mazer (Librairie polytechnique Ch. 
Béranger, Paris 1955). 280 pp., 159 fig.; relié fr. fr. 4975,—. 

De nombreux ouvrages ont été publiés sur la Mécanique vibratoire. Celui de 
M. Maze, professeur a la Faculté des Sciences de Poitiers, reproduit les lecons 
données par cet auteur a l’Ecole nationale supérieure de l’Aéronautique, de 
1950 a 1954. ' 

Il s’agit avant tout de l’étude des systémes mécaniques dits «linéaires», dont 
lévolution reléve, avec une précision suffisante, d’équations différentielles 
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linéaires a coefficients constants. Seul le dernier des huit chapitres dont voici tes 
titres fait exception: 

Généralités sur les systémes vibrants. — Systémes vibrants conservatifs a un 
degré de liberté. — Systémes vibrants dissipatifs A un degré de liberté. — Notion 
de couplage. Systémes conservatifs &4-deux degrés de liberté et plus. — Systémes 
en chaine. Filtres. Systémes continus conservatifs. — Systémes dissipatifs a 
plusieurs degrés de liberté. — Systémes vibrants associés A un réservoir d’énergie. — 
Intervention des facteurs non linéaires. 

M. Mazer a cherché, dans ces legons, a atteindre deux buts: Tout d’abord, 
exposer la dynamique des systémes linéaires, en mentionnant constamment les 
analogies entre ces systémes et les phénoménes appartenant a d’autres branches 
de la physique (€lectricité, acoustique, etc.). Ensuite, traiter le plus grand nombre 
possible d’exemples intéressant la technique aéronautique. Il y a d’autant mieux 
réussi que ses fonctions de directeur scientifique de l’Office national d’Etudes et 
de Recherches aéronautiques lui ont permis d’acquérir une grande expérience: 
dans ce domaine. 

Nous pouvons recommander vivement la lecture de cet ouvrage et nous ne 
doutons pas qu’il «rende de bons et loyaux services, non seulement a des étudiants 
ou a des ingénieurs débutants, mais aussi, et non moins, a des ingénieurs plus 
avancés», comme le dit dans la préface M. MaurIcE Roy, membre de l’Académie 
des Sciences. H. Favre 


Praxis der Differentialgleichungen. Von H. v. SANDEN, 4. Aufl. (W. de 
Gruyter, Berlin 1955). 114 S., 21 Abb.; DM. 6.80. 

Die vorliegende Schrift, aus Vorlesungen des Verfassers entstanden, liefert 
auf gedrangtem Raum die Grundlage zur numerischen Integration gewohnlicher 
Differentialgleichungen. Der Verfasser beschrankt sich auf Differentialgleichun- 
gen 2. Ordnung (und Systeme von solchen) sowie auf die wichtigsten Methoden 
(Apams und RunceE-Kutta). Weiter werden auch Rand- und Eigenwertauf- 
gaben behandelt, wobei unter anderem auch das Ritzsche Verfahren zur Sprache 
kommt. Das kleine Werk kann gerade, weil es sich auf das Wesentliche beschrankt, 
als Lehrbuch zur Einfiihrung in die Materie bestens empfohlen werden. 

H. Rutishauser 


Transactions of Symposia in Applied Mathematics. Band 2 (Inter- 
science Publishers Inc., New York 1955). Communications on Pure and Applied 
Mathematics 8, 1 (Februar 1955; Sonderheft). $5.00. 

Das Buch vereinigt die an einer Tagung an der University of Chicago am 
29./30. April 1954 gehaltenen Vortrage: 1. P. M. Morse, Operations Research. — 
2. Jerzy Newman, The Problem of Inductive Inference. — 3. H. O. HARTLEY, 
Some Recent Developments in Analysis of Variance. — 4. J. E. Maver, Two Un- 
solved Problems in Statistical Mechanics. — 5. M. R. HestEenes, Iterative Compu- 
tational Methods. — 6. JoHN Topp, Motivations for Working in Numerical Analysis. 
— 7. A.A. BENNETT, Some Numerical Computations in Ordnance Problems. — 
8. C. A. TRUESDELL, The Simplest Rate Theory of Pure Elasticity. — 9. J. J. Sto- 
KER, On the Stability of Mechanical Systems. — 10. FLORENT BUREAU, Divergent 
Integrals and Partial Differential Equations. — 11. WILLIAM FELLER, On Differential 
Operators and Boundary Conditions. ie 

Die behandelten Themen geben einen recht guten Uberblick tiber die gegen- 
-wartig im Vordergrund des Interesses stehenden Probleme der angewandten 
Mathematik. HA. Rutishauser 
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Scattering and Diffraction of Radio Waves. Von J. K. MENTZER (Per- 
gamon Press Ltd., London 1955). 134 S., 34 Fig.; £1.10.0. : 

Die wissenschaftliche Biicherfolge der Pergamon Press, «Electronics and 
Wave Series», herausgegeben von D. W. Fry, hat in der Fachwelt rasch Beach- 
tung gefunden. Der 7. Band, Scattering and Diffraction of Radio Waves, reiht 
sich den friiheren in gliicklicher Weise an. Die Streuungs- und Beugungsprobleme 
der Radiowellen, welche namentlich im Gebiete kurzer und kitirzester Wellen 
sehr zahlreich sind, bieten bekanntlich mancherlei Schwierigkeiten, was fiir die 
theoretische Erfassung vieler Fragen von grossem Nachteil ist. Eine sorgfaltig 
begriindete und systematische Zusammenstellung der Methoden nebst Behand- 
lung der typischen Beispiele verschiedener Art wird daher weiten Fachkreisen 
willkommen sein. Diese Aufgabe ist durch MENTZER in seinem Buche sehr gut 
gelést worden. Die Darstellung ist tibersichtlich und klar und durch zahlreiche 
Figuren erganzt. Es liegt in der Natur der, Sache, dass ein gewisses Mass von 
Mathematik nicht vermieden werden kann und dass das Buch studiert sein will. 
Aber das Studium lohnt sich; wir empfehlen es jedem, welcher auf diesem Fach- 
gebiet wissenschaftlich tatig ist. F. Tank 


Hochvakuum-Elektronenréhren. Band 1: Physikalische Grundlagen. 
Von H. RorHE und W. KLEEN (Akademische Verlagsgesellschaft m. b. H., Frank- 
furt am Main 1955). 297 S., 182 Abb.; DM 32.-. 

Das vierbandige Werk von H. RotHEe und W. KLEEN iiber die Grundlagen 
und Anwendungen der Elektronenr6éhren ist in weiten Kreisen wohlbekannt. 
Da dessen Entstehung iiber 15 Jahre zuriickliegt, kann eine Umarbeitung und 
Anpassung desselben an die Fortschritte der neuesten Zeit nur sehr willkommen 
geheissen werden. Der erste Band, Hochvakuum-Elektronenrohren, Physikalische 
Grundlagen, liegt nunmehr in neuer Fassung vor. Der reiche Stoff ist sorgfaltig 
ausgewahlt und in itibersichtlicher, fasslicher Weise dargestellt. Insbesondere 
ist auch den so wichtigen Teilgebieten, wie Kenntnis der Potentialfelder mit und 
ohne Raumladung, Ermittlung von Elektronenbahnen (Elektronenoptik), Elek- 
tronenlaufzeiten usw., eine eingehende Berticksichtigung geschenkt. Die ausfiihr- 
lichen Literaturangaben am Schlusse eines jeden Kapitels werden vielen Lesern 
sehr willkommen sein. Druck und Figuren sind von guter Qualitat. Das Werk 
sei in seiner neuen Gestalt Studierenden, Physikern und Ingenieuren, welche 
sich ein vertieftes Wissen tiber die Physik der Hochvakuum-Elektronenréhren 
aneignen wollen, auf das beste empfohlen. F. Tank 


Schaltungstheorie und Messtechnik des Dezimeter- und Zentimeter- 
wellengebietes. Von A. WeEIssFLocH (Birkhauser Verlag, Basel und Stuttgart 
1954). 308 S., 288 Abb.; geb. Fr./DM 33.50, brosch. Fr./DM 29.30. - 

Die Erscheinungen im Gebiete der elektrischen Dezimeter- und Zentimeter- 
wellen besitzen bekanntlich ausgesprochenen Wellencharakter, weshalb man zu 
ihrer Behandlung auf die Maxwellschen Gleichungen und deren Lésungen bei 
bestimmten Randbedingungen (Hohlleiter, Hohlraumresonatoren usw.) zuriick- 
zugreifen pflegt. Sind die Grundphanomene einmal erfasst, so zeigt es sich, dass 
in der Schaltungs- und Messtechnik linear gebrochene Beziehungen (Transfor- 
mation, Anpassung usw.) eine wichtige Rolle spielen, in formal durchaus ahnlicher 
Weise wie in der Theorie der Leitungen und der linearen Netzwerke (Vierpoltheorie 
usw.). Dem Vorwort entnehmen wir die folgenden Satze: «Die vorliegende Dar- 
stellung macht keinen Gebrauch von den Maxwellschen Gleichungen. Zur exakten 
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Berechnung, zum Beispiel der Feldverteilung in Hohlleitern, wird auf die zahlreichen 
existierenden Werke verwiesen... Als mathematisches Werkzeug tritt jedoch 
die Geometrie der linear gebrochenen Funktionen stark in den Vordergrund. Es 
ist dies ein mathematischer Zweig mit sehr einfachen und anschaulichen Ge- 
setzen...» Wer sich mit diesen niitzlichen, vielseitig anwendbaren und allgemei- 
nen theoretischen Betrachtungen vertraut machen will, dem sei das Studium des 
Buches bestens empfohlen. Die Darstellung ist ausfiihrlich und klar, die Ausstat- 
tung des Buches durch den Verlag vorbildlich. F. Tank 


An Introduction to Reactor Physics. Von D. J.Lirrter und J.F. RAFFLE 
(Pergamon Press Ltd., London 1955). 196 S., 80 Fig.; 25s. 

Die knappe Einfiihrung ist aus den Vorlesungen, welche der Reaktorschule 
in Harwell zugrunde liegen, hervorgegangen. Sie wendet sich an einen Leser- 
kreis ohne Vorkenntnisse in Kernphysik. Da Harwell vor allem den Bau und 
Betrieb von Reaktoren mit Graphitmoderatoren geférdert hat, sind die im vor- 
liegenden Bandchen als Erlauterung verwendeten Beispiele auf diese Art von 
Reaktoren abgestimmt. Etwa ein Viertel des Inhalts dient der Vermittlung der 
fiir das Verstandnis der Reaktorphysik unumganglichen kernphysikalischen 
Grundbegriffe. Der reaktortheoretische Teil bedient sich einer von der ameri- 
kanischen Literatur abweichenden Schreibweise, indem in den Formeln an Stelle 
des Neutronenflusses (Neutronen pro Quadratzentimeter und Sekunde) iiberall 
die Neutronendichte (Neutronen pro Kubikzentimeter) verwendet wird. Konse- 
quenterweise erscheint dann auch die in der Diffusionstheorie der Gase gebrauch- 
liche Diffusionskonstante der Dichte. Diese Schreibweise erschwert dem Anfanger, 
der gleichzeitig zu amerikanischen Verdffentlichungen greift, das Studium. 

Das letzte Viertel des Buches enthalt Hinweise tiber Aufbau und Arbeitsweise 
von im Reaktorbetrieb gebrauchlichen Strahlendetektoren sowie einen Abschnitt 
tiber Effekte, welche im Kristallgefiige eines Festkérpers durch Bestrahlung auf- 
treten. 

Ein Anhang vermittelt die heute gebrauchlichen Definitionen und Toleranz- 
werte fiir die Einwirkung von Strahlungen auf den lebenden Organismus. 

W. Haig 


Praktische Mathematik. Von H. von SanpEN (B. G. Teubner Verlags- 
gesellschaft, Stuttgart 1956 [Teubners Mathematische Leitfaden, Band 44)). 
Toro, SO ADD. DM: 760. 

Das friiher hier schon einmal besprochene Bandchen aus der Reihe der Teub- 
nerschen mathematischen Leitfaden ist jetzt in vierter, iiberarbeiteter Auflage 
herausgekommen. Den meisten Kapiteln wurden Erganzungen beigefiigt, so iiber 
die Approximation der Ableitung einer Funktion, quadratische Interpolation 
nicht 4quidistanter Funktionswerte, das Verfahren von GRAEFFE, ferner einige 
Beispiele zur Nomographie. Im iibrigen wurde der bisherige Aufbau beibehalten: 
Darstellung von Funktionen; Satz von Taytor; Auflosung von Gleichungen ; 
Integration, Differentiation und Interpolation; Statistik; Ausgleichsrechnung : 
harmonische Analyse und Synthese; Nomographie. Durch diese Erweiterungen 
diirfte diese fiir Studenten der Ingenieurwissenschaften und ebenso fiir den 
Praktiker gedachte Einfiihrung an Brauchbarkeit noch gewonnen haben. — Nach- 
zutragen ist noch, dass — wie aus uns vorgelegten Unterlagen hervorgeht = der 
- neue Sitz des Teubnerschen Verlages sich in Stuttgart und nicht mehr in Leipzig 
befindet. E. Roth-Desmeules 
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Vorlesungen iiber Orthogonalreihen. Von F. Tricomr (Springer-Verlag, 
Berlin 1955). 264 S., 13 Abb.; DM 34.—/37.60. 
Aus der Theorie der Fourier-Reihen ist als Verallgemeinerung die Theorie der 
Orthogonalfunktionen, das heisst der Funktionssysteme 9,(¥) (vy = 0, 1, 2, ...) 
b 


mit der Eigenschaft ’ P(*) pul) P(~) dx = 0 (uw + 1) herausgewachsen, Diese 


a 
Theorie hat als wichtiges Hilfsmittel der mathematischen Physik in den letzten 
Jahren noch stark an Bedeutung gewonnen. Das Erscheinen eines Lehrbuches 
iiber diesen Gegenstand ist daher warmstens zu begriissen; um so mehr, als eine 
geschlossene Darstellung der Theorie in der deutschsprachigen Literatur bisher 
eigentlich gefehlt hat. 

Inhalt. Kapitel I: Allgemeines tiber orthogonale Funktionssysteme — II: 
Theorie der trigonometrischen Reihen — III: Konvergenzeigenschaften der trigo- 
nometrischen Reihen — IV: Allgemeines iiber orthogonale Polynome. V, VI: Klas- 
sische orthogonale Polynome mit endlichem bzw. unendlichem Grundintervall. 

H. Rutishauser 


Selecta Hermann Wey] (Birkhauser Verlag, Basel und Stuttgart 1956). 592S.; 
Fr. (DM) 48.90. 

Am 8. Dezember 1955 ist HERMANN WEvL gestorben. Wenige Wochen vorher, 
an seinem 70. Geburtstag, wurde ihm zu seiner grossen Freude das erste Exemplar 
seiner Selecta iiberreicht, bei deren Auswahl er selber mitgewirkt hat. So sind 
denn die schénsten und wichtigsten Arbeiten WEYLS zu einem stattlichen Bande 
vereinigt. Alle sind Kunstwerke, der Sprache und dem Inhalte nach. 

An der Spitze steht die beriihmte Abhandlung von 1910 tiber gewohnliche 
Differentialgleichungen mit Singularitaten, in der zum ersten Male prdazisiert 
wird, wie man die Randbedingungen des Eigenwertproblems einer gewohnlichen 
Differentialgleichung 2. Ordnung anzusetzen hat, damit willkiirliche Funktionen 
sich nach Eigenfunktionen entwickeln lassen. An diese klassische Abhandlung 
schliessen sich weitere, ebenfalls grundlegende Arbeiten tiber Eigenwertprobleme 
der Physik und Mathematik. Darunter fallen die Arbeit von 1925 iiber Integral- 
gleichungen und fastperiodische Funktionen, die beriihmte dreiteilige Abhand- 
lung von 1925-26 tiber Darstellung der halbeinfachen Lieschen Gruppen und die 
daran anschliessende gemeinsame Arbeit mit F. PETER iiber geschlossene Liesche 
Gruppen. 

Die Grundlagendiskussion und die Arbeiten iiber Algebra und Zahlentheorie, 
Differentialgeometrie, Funktionentheorie und harmonische Integrale muss ich 
leider hier beiseite lassen, ebenso die iiber Gravitation und Elektrizitat und iiber 
das Symmetrieproblem der Quantenmechanik, um mich auf die Beitrage zu be- 
schranken, welche fiir die angewandte Mathematik und Physik von Wichtigkeit 
sind. 
Da ist zunachst die Abhandlung von 1915 iiber die asymptotische Verteilung 
der Eigenschwingungen eines elastischen Kérpers. Sodann die sehr schéne Arbeit 
von 1919 iiber die Ausbreitung elektromagnetischer Wellen iiber einem ebenen 
Leiter, in der ein wichtiges Problem der Radiotelegraphie, das SoMMERFELD 
zuerst behandelt hatte, vollstandig gelést wird. An SoMMERFELDS Untersuchun- 
gen kniipft auch die letzte Arbeit iiber die natiirlichen Randwertaufgaben im 
Aussenraum fiir Strahlenfelder an. pe 

Sehr anregend ist die Abhandlung von 1940 iiber die Methode der orthogo- 
nalen Projektion in der Potentialtheorie. 
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Praktische Fragen der Hydrodynamik bilden den Ausgangspunkt der Arbeit 
von 1942 tiber die Differentialgleichungen von Grenzschichtproblemen. Die-Dif- 
ferentialgleichungen werden durch gut konvergente Reihen explizit gelost. 

In der immer anschwellenden Flut der Fachliteratur ist es sehr wichtig, die 
wertvollsten Arbeiten der wahrhaft grossen Forscher in Buchform vereinigt zu 
haben. Zu diesen Grossen gehérte WEYL. B.L. van der Waerden 


Frequency Response. Herausgegeben von R. OLDENBURGER (Macmillan 
Company, New York 1956). 392 S., 404 Abb.; $7.50. 

In diesem Buche werden vor allem die an der Tagung Frequency Response 
Symposium vom Dezember 1953 in New York prasentierten Aufsaitze wieder- 
gegeben, gemeinsam mit weiteren Arbeiten im Zusammenhang mit der Frequenz- 
gangmethode. 

Die 28 von verschiedenen Verfassern stammenden Aufsdtze beschlagen die 
unterschiedlichsten Gebiete der Theorie und Anwendung des Frequenzgangver- 
fahrens. Eine erste Gruppe von Aufsatzen ist den allgemeinen Grundlagen sowie 
der — was europaische Quellen anbelangt, allerdings liickenhaften — Bibliographie 
gewidmet. Weitere Arbeiten betreffen die experimentelle Bestimmung des Fre- 
quenzganges und die dazu bendétigten Hilfsmittel, namentlich Schwingungser- 
zeuger. Zwei Gruppen von Aufsadtzen gelten der Anwendung auf Probleme der 
Praxis. Unter anderem werden Beispiele aus den Gebieten der Regelung von 
Turbinen, Flugzeugen und chemischen Prozessen behandelt. Andere Arbeiten 
gelten vorwiegend der Weiterentwicklung der Frequenzgangmethodik und ihrer 
Anpassung an besondere Problemstellungen der Praxis. So betreffen zwei Grup- 
pen von Aufsadtzen den Zusammenhang zwischen Frequenzgang und Ubergangs- 
funktion bzw. Verfahren zur Optimierung von Reglerkennwerten. Zwei weitere 
Gruppen sind der Ubertragung der Frequenzgangtheorie auf nicht lineare bzw. 
unstetige Systeme gewidmet. Schliesslich behandelt auch ein Aufsatz die An- 
wendung statistischer Methoden in diesem Zusammenhang. 

Das Buch vermittelt in seiner Vielseitigkeit eine Fiille von theoretischen und 
praktischen Erkenntnissen und wird vor allem dem Regelungsfachmann ausge- 
zeichnete Dienste leisten. Manchem nichtamerikanischen Leser mag die durch die 
grosse Zahl der Autoren bedingte Uneinheitlichkeit der Begriffsbildung und Dar- 
stellungsweise die Lektiire etwas erschweren. Anderseits wird er gerade durch 
die damit zusammenhangende Vielfalt und Farbigkeit der Aspekte reichlich 
entschadigt werden. P. Profos 


An Introduction to the Theory of Aeroelasticity. Von Y. C. Fune (John 
Wiley & Sons, New York 1955). 490 S., 143 Fig.; $10.50. 

Das Problem der aerodynamisch angeregten Schwingungen ist sicher eines 
der wichtigsten und schwierigsten des Flugzeugbaues. Aber auch in anderen 
Gebieten der Technik sind diese Fragen aktuell, so dass eine gute Einfiihrung in 
die Aeroelastizitaét zweifellos einem Bediirfnis entspricht. Das vorliegende Buch 
bringt eine klare, logisch aufgebaute Darstellung des sehr ausgedehnten Frage- 
komplexes, mit spezieller Erlauterung der physikalischen Grundlagen des 
Schwingungsvorganges. Die verwendeten mathematischen Methoden werden von 

Fall zu Fall knapp, aber zum Verstandnis geniigend genau behandelt. 

Die ersten Kapitel sind der Flugzeugstatik, der Aerodynamik und allgemeinen 
Schwingungs- und Instabilitatserscheinungen gewidmet, als Vorbereitung zur 
‘Behandlung der eigentlichen Fliigelschwingungen mit mehreren Freiheitsgraden. 
Der Einfluss der atmosphiarischen Turbulenz und die AbreiSschwingungen sind 


a 
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ebenfalls Gegenstand einer eingehenden Diskussion. Von besonderem Wert ist die: 
ausgedehnte Bibliographie. Das Buch ist sehr fliissig geschrieben; die mathema- 
tischen Ableitungen sind knapp, aber streng gehalten, die numerischen Beispiele 
sind vielleicht zu wenig zahireich, um das notwendige Gefiihl fiir die Gréssenord- 
nungen zu erwecken. 

Im grossen ganzen ein ausgezeichnetes Werk, das den studierenden Flug- 
ingenieuren bestens empfohlen werden kann. P. de Haller 


Hédhere Mathematik, Teil VII: Raumliche und ebene Potentialfunktionen. 
Konforme Abbildung. Integralgleichungen. Variationsrechnung. Von R. RoTHE, 
herausgegeben von W. SCHMEIDLER (B. G. Teubner Verlagsgesellschaft, Stutt- 
gart 1956). 218 S., 43 Abb.; DM 19.80/Fr. 22.40. 

Dem von Rupotr Rorue begriindeten und nun von W. SCHMEIDLER heraus- 
gegebenen Unterrichtswerk der héhern Mathematik fiir Mathematiker, Physiker — 
und Ingenieure wurde jetzt noch ein siebenter und abschliessender Band hinzu- ~ 
gefiigt, der im wesentlichen drei Gebiete umfasst: Potentialfunktionen, Integral- — 
gleichungen und Variationsrechnung. Der erste Teil behandelt die raumlichen 
und ebenen Potentialfunktionen und bringt erganzend noch einiges aus der kon-_ 
formen Abbildung (Polygonabbildung) und der Funktionentheorie (elliptische 
Funktionen). Dann folgen im zweiten Teil eine determinantenfreie Darstellung 
der linearen Gleichungssysteme und Matrizen als Vorbereitung fiir die Lésung 
der Integralgleichungen erster und zweiter Art. Schliesslich sind im letzten Teil 
noch einige Ausfiihrungen zur Variationsrechnung enthalten, insbesondere tiber — 
Probleme aus der Differentialgeometrie der Flachen und die WeierstraBsche 
sowie die Jacobi-Hamiltonsche Theorie. — Naturgemass zwang der vorhandene — 
Raum zu einer betrachtlichen Stoffbegrenzung und entsprechend zu einer in 
geschickter Weise getroffenen Auswahl der zu behandelnden Fragen. | 

Wie bei den iibrigen Banden wurde auch diesmal jedem Abschnitt eine — 
Reihe Aufgaben (durchschnittlich sind es deren 20) beigefiigt; fermer sind Hin-~ 
weise auf ausfiihrlichere Werke vorhanden. Es ist nicht zu zweifeln, dass diese . 
nun abgeschlossene, den heutigen Erfordernissen weitgehend geniigende Ein- — 
fiihrung in die hédhere Mathematik und ihre Anwendungen vor allem ftir die 
Studierenden der Mathematik und der Physik eine wichtige und zuverlassige — 
Hilfe sein wird, zumal die die Erfahrungen des Verfassers widerspiegelnde Dar- | 
stellung sehr sorgfaltig und iibersichtlich ist. E. Roth-Desmeules 


: 
: 


The Constitutional Diagrams of Alloys. A Bibliography. Urspriinglich — 
zusammengestellt von J. L. Haucuton, D. Sc. Zweite Auflage zusammengestellt 
von A. PRINCE und B. Mer., A. I. M. (The Institute of Metals, London 1956). 
3235. sEeG1 

Die erste Ausgabe dieser in ihrer Vollstandigkeit wohl einzigartigen Biblio-— 
graphie erschien 1942. Seither sind die Phasengrenzen zahlreicher weiterer Zwei-, _ 
Drei- und Vierstofflegierungen erforscht und veréffentlicht oder friihere Erkennt- 
nisse verbessert worden. Die 2. Auflage dieses fiir den Metallkundler unentbehr- 
lichen Hilfsbuches enthalt schatzungsweise 12000 bis 13000 Quellenangaben iiber 
diese Arbeiten in allen Landern, die bis Mitte 1955 nachgefiihrt sind. Die Ubersicht 
ist durch die streng alphabetische Ordnung der Systeme an Hand der Symbole | 
fiir die chemischen Elemente sehr erleichtert. E. Bickel 
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I. Das Wesen des Transformators ~ II. Die Gestalt des Transformators — III. Das Wesen des Drehstrom- 

transformators — IV. Die Warmeabfuhr aus dem arbeitenden Transformator — V. Der Eisenkern. Der 

Leerlauf. Der Einschaltstromstoss = VI. Die Wicklungen — VII. Randprobleme des Transformatorenbaues 
VIII. Der Entwurf. Beispiele von ausgefiihrten Transformatoren. 


Die Transformatorenprobleme haben eine verhdltnismdssig gewaltige Vergangenheit hinter sich, sind 
indessen immer noch nicht vollstandig gelést. Es gibt eine Unmenge von mitwirkenden Faktoren, die den 
Bau eines grossen Transformators beeinflussen, und hinter jedem Faktor steht ein Problem. 

«Die Transformatoren», deren dritte Auflage nun herauskommt, haben durch den Erfolg ihrer beiden 
vorangegangenen Auflagen bewiesen, dass sie die Problematik des Transformatorenbaues richtig ange- 
packt und zufriedenstellend behandelt haben. Seit dem Erscheinen der zweiten Auflage sind allerdings 
mehr als 30 Jahre vergangen, so dass die dritte Auflage den Rahmen der zu behandelnden Probleme 
ganz erheblich erweitern musste. Sie hat vor allem die Aufgaben zu behandeln gehabt, die mit den sehr 
hohen Spannungen im Transformatorenbau aufgetreten sind. Dariiber hinaus gab es neue konstruktive 
Probleme, die zum Teil mit den ganz grossen Leistungen zusammenhangen, zum Teil mit dem Auftauchen 
neuer Baustoffe verbunden sind. 

Dies alles soll die dritte Auflage «Der Iransformatoren» dem Sonderfach, in das sie ausgerichtet ist, 
bringen. Sie kommt trotz der erweiterten Problematik mit einem etwas bescheideneren Umfang aus als die 
vorangegangene zweite Auflage. Trotzdem bemiiht sie sich, die praktische Anwendung der beschriebenen 
Problemlésungen anhand von Beispielen zu erleichtern. 
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